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ABSTRACT 

Using a simple description of feedback from black hole growth, we develop an analytic model for the fueling 
of Seyferts (low-luminosity AGN) and their relation to their host galaxies, Eddington ratio distributions, and 
cosmological evolution. We derive a solution for the time evolution of accretion rates in a feedback-driven 
blast wave, applicable to large-scale outflows from bright quasars in galaxy mergers, low-luminosity AGN, and 
black holes or neutron stars in supernova remnants. Under the assumption that cold gas stochastically accretes 
onto a central supermassive black hole at a rate set by the dynamics of that gas, our solution determines the 
evolution of Seyfert light curves. Using this model, we predict the Seyfert luminosity function, duty cycles and 
AGN "lifetimes," and the distribution of host morphologies, Eddington ratios, and obscuration as a function 
of AGN luminosity and black hole mass, and find agreement with observations at z = 0. We consider the 
breakdown of the contribution from this mechanism and from stellar wind and virialized hot gas accretion and 
merger-driven activity. We also make specific predictions for the weak evolution of the Seyfert luminosity 
function; i.e. luminosity function of quiescent, as opposed to merger-driven activity, as a function of redshift, 
and for changes in both the slope and scatter of the M^n-cr relation at Iow-Mbh- Our modeling provides a 
quantitative and physical distinction between local, low-luminosity quiescent AGN activity and violent, merger- 
driven bright quasars. In our picture, the quiescent mode of fueling dominates over a wide range of luminosities 
(-14 > Mb > -22) at z = 0, where most black hole growth occurs in objects with Mbh ^ 10^ Mq, in SO and 
Sa/b galaxies. However, quasar activity from gas-rich mergers evolves more rapidly with redshift, and by z = 1, 
quiescent fueling is important only at luminosities an order of magnitude or more below the "break" in the 
luminosity function. Consequently, although non-merger driven fueling is important for black hole growth 
and the Mbh - o" relation at low Mbh, it does not significantly contribute to the black hole mass density of the 
Universe or to cosmological backgrounds. 

Subject headings: quasars: general — galaxies: active — galaxies: evolution — cosmology: theory 
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1. INTRODUCTION 

While it is now generally accepted that quasars and 
active galactic nuclei (AGN) are powered by the accre- 
tion of gas onto supermassive black holes in the centers 
of galaxies (e.g. Salpeter 1964; Zel'dovich & Novikov 
1964; Lynden-Bell 1969), the mechanism that fuels these 
objects is uncertain. In the local Universe, ultralu- 
minous infrared galaxies (ULIRGs) have bolometric lu- 
minosities similar to bright quasars a nd ar e always in 
mergers (see e.g., 'Sanders et all 119861 Il988t ISoifer et al.1 
[1987; Sanders & Mirabel 199^rand multi-wavelength sfiid- 
ies have shown that man y contain growing , optically 
obscured black holes (Komossa.et_alJ 2003; Gerssen et al 



|2p04; Max et al. 2005; Alexander et al. 2005a b; Bor vs et al 
|2005). Observations of low-redshift quasar hosts also 
reveal a connection bet ween galaxy mergers and quasar 
activity (e.g., Heckman et al. 1984; Stockton & RidswaV 
[1991; Hutchings & Neff 1992; BahcaU et al. 1996 , 1997 ; 
[Canalizo & Stocktoni lioOll) . These facts motivate a scenario 
where mergers of gas-rich galaxies provide the fuel to power 
nuclear starbursts which evolve into bright quasars. 

Numerical simulations have identified a physical mecha- 
nism for transporting gas into the inner regions of galax- 
ies in a merger by showing that tidally-induced gravitational 
torques remove angular momentum from the gas (e.g. Barnes 
& Hernquist 1991, 1996; Mihos & Hernquist 1996). Models 

' We dedicate this paper to lohn Balicall, who inspired our interest in this 
problem. 

- Hai'vard-Smithsonian Center for Astrophysics, 60 Garden Street, Cam- 
bridge, MA 02138. USA 



including supermassi ve black holes (iDi Matteo et al. 2005[ 
ISpringel et alJ l2005b() support the conjecture that ULIRGs 
evolve into quasars and suggest that feedback from black hole 
growth mediates this transition by expelling obscuring gas 
and dust. This " blowout" results in a shor t-lived, bright opti- 
cal quasar ( Hop kins etalJ2005albt l2()06i and eventually ter- 
minates the activity, leaving a remnant that quickly reddens 
(e.g. Springel et al. 2005a) and satisfies obs erved correla- 
tions between black hole mass and the mass ( Magorria n et alJ 
[1998; McLure & Dunlop 2002; Marconi & Hun t 2003) or ve- 
locity dispersion (i.e. the Mbh-c relation: Ferrarese & Merritt 
2000; Gebhardt et al. 2000) of spheroids. 

Moreover, the time evolution in these simulations repro- 
duces many quasar observables, including luminosity func- 
tions, host prop erties, and accretion rates ( Hopkins et all 
l2005cldil2()06aB . The modeling also indicates that the bright 
quasar population and the bulk of the cosmic quasar lumi- 
nosity density and buildup of the black hole mass d ensity 
must be dominated by merger-driven growth ( Hopkins et ajj 
|2005e), to be consistent with the evolution of the X-ray back- 
ground ( Hopkins et al. 20 06al) . th e red galaxy age, metal licity, 
mass, and luminosity functions JHopkins et alJl2006lA . and 
merge r luminosity fun ctions and star formation rate densities 
( Hopk ins et al.l2006el) . 

However, many local, low-luminosity AGN resi de in 
quiesc e nt, non-interacting galaxies ( e.g., Kauffmann^^^ 
l2003at 'Sanchez & Gonzalez-Serran ol 12003; .Sanchez et alJ 
l2004t iMcLure & Dunlop 2004; Ha o et al.ll2005h . Some of 
these objects do in fact fit naturally into a merger-driven pic- 
ture. At relatively high luminosities, ellipticals with young 
stellar populations and moderate Eddington ratios have ob- 
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served properties consistent with decay in the quasar light 
curve from a pr evious bright quasar epoch in a spheroid- 
forming merger ("Kauffmann et alJ l2003at iBest et all l2005t 
[Hopkins et al. 2005e, 2006c). At the lowest luminosities char- 
acteristic of low-luminosity AGN and LINERs, populations 
of very low accretion rate "dead" ellipticals dominate, fu- 
eled via accretion of hot (viri alized) spheroid gas and steady 
mass loss from sta rs (e.g., ICiotti'& O strikei^ '1997", '20011 
l^llegrini 2005; So ria et alJ l2005bl) . presumably in radia- 
tively inefficient accretion states (e.g., Naravan & Yilll995h . 
Although the detailed prope rties of the se objects are sub- 
ject to debate (see, e.g., NaravanI 120041) . the fuel source is 
reasonably well-understood, and it is clear from compar- 
ison o f black hole mass and quasar luminosity functions 
(e.g., Soltan 1982; Salucci et al. 1999; Yu & Tremaine 200l 
iMar coni et al. 2004; Shankaretal. 2004), background syn- 
thesis models (e .g., Comastri et al.l |1995l: iGilli etalJ 119991 
lElvis et alJ [20021 "iJeda et al. 2003; "Cao"2005'), and obser- 
vations of the accretion rate distribution (Vestergaard 2004j^ 
iHeckman et al. 2004; Mc Lure & Du nlop 2004; Hopkin s et all 
|2p05e) that these modes do not dominate black hole growth 
or cosmological backgrounds. 

Nevertheless, a number of relatively high accretion rate 
objects are observed at low reds hift in undisturbed, late- 
type, star-forming galaxies (e.g., iKauffma rm et al. 2003d) 
with small (Mbh < 10^ M©) black holes (HeckmanetaQ 
|2004). Indeed, the local AGN population largely comprises 
black holes in non-interacting, star-forming SO and Sa/b hosts 
with n o evidence of galaxy-scale perturbations or distur- 
bances ("Dong & De Robertis'2005"), spanning most of the ob- 
served z = Q AGN luminosity function (Hao et al. 2005). Al- 
though these AGN have been more thoroughly studied and 
are more well-understood than bright quasars at high red- 
shift, there is no self-consistent model for their triggering, 
fueling, and evolution. Previous investigations of their fu- 
eling have mainly been restricted to estimating whether or 
not a giv en fuel source could provide an adequat e mass sup- 
ply (e.g.. ' Lvnden-Bell.1969: .Hills.. 1975.; .Shields & Wheeled 
[1978; Mat hewslll99'ol) . or have examined the evolution of 
light curves in limited subclasses of these objects (e.g., 
[Ciotti & Ostriker 2001) . A critical uncertainty in previous ef- 
forts to model AGN is the light curve evolution once a given 
"trigger" occurs. 

Independent of how fuel is delivered to the supermassive 
black holes in Seyferts and low-luminosity AGN, their subse- 
quent evolution may resemble that in bright quasars if feed- 
back from accretion regulates black hole growth. In the sim- 
ulations of Di Matteo et al. (2005), the impact of this feed- 
back resembles an explosion because the energy is deposited 
on scales small compared to the host galaxy and because the 
black hole grows nearly exponentially during the quasar phase 
with an e-folding time that is short compared to the character- 
istic dynamical time of the host potential. Indeed, Hopkins et 
al. (2005g) show that the outflow driven by this feedback can 
be described as a blast wave. In what follows, we investigate 
whether solutions of this type can also be used to characterize 
the evolution of low-luminosity AGN and Seyferts. 

There is observation evidence for outflows and winds in 
both AGN and quasars (for a review, see Veilleux et al. 2005jj). 
The kinematics of gas in the narrow line regions of local, 
low-luminosity AGN (see, e.g. Rice et al. 2006, and refer- 
ences therein) display bi-conical or nearly isotropic (wide- 
angle) radial outflow s at speeds ^ 10^ - lO-'kms"' (e.g., 
iCrenshaw et alJl2000l) . with wind dynamical times ~ 10 - 



10^ y r and entrained molecu lar gas mas ses ~ 10^ - 10''Mo 
(e.g.. iStark & CarlsoiJlT984t IWalter et al] |2002l) . The "warm 
absorber" seen in soft X-rays is also indicative of a significant 
outflow generated local to the AGN (e.g., Laor et al. 1997, and 
references therein), accelerating radially to larger terminal 
veloc ities > afewlOOkms"' (Ruiz et al, 2001; Kas pi et alJ 
l2002h . These absorbing outflows appear ubiquitous (and 
are likely even where not directly observed, given ob- 
served wind covering angles and clumping factors; see e.g. 
iRupke et al. ' '2005 ) and ar e associated with clumpy, high- 
ionization structures (e.g., ICrenshaw et al . 1999). Outflow 
energe tics and entrained masses appear to scale with AGN 
power (iBaum & McCa rthv'20001). and intense winds with ve- 
locities ~ lO'^kms"' are seen i n bright, broad absorption line 
(BAL) quasars (e.g.. Reeves et al.'2003"). Observations of typ- 
ical, narrow-line quasars also find BAL-like outflow velocities 
for som e species acce lerating at small radii from the central 
engine (iPounds et al.ll2003a.b.) . and outflows with large cov- 
ering angles in central regions of "normal" quasars may eve n 
be detected in gravitational lensing signatures (Green 2006). 

Although the relative contributions of star formation and 
AGN to the energetics i s unclear (see e.g., Baum et al] 
ll993tlLevenson et alJ2001l) . high-resolution data indicate that 
even Seyfert II ULIRGs in which large-scale winds may be 
driven by star formation have central, AGN-driven outflows 
([Cecil et al. 2001; Rupke et al. 2005), and observations at op- 
tical. X-ray, and radio w avelengths (Colbert et al. 1996a bj 
119981 IWhittle & Wilsonl[2004<) identify a number of mech- 
anisms by which AGN power can couple to outflows over 
'-^kpc scales. The decomposition of the narrow line region 
components also suggests wind ram pressure or radiation 
pressure from the central source as a driver dKaiser et alJ 
.2000i) . and in some cases the wind injection zone is ob- 
served t o be small relative to the scale of nuclear star for- 
mation ("Smith & Wilson"2001'). Furthermore, similar out- 
flow structures are observed in the hot gas of some spheroidal 
systems wh ich have no rapid associated star formation (e.g., 
iBiller et aU2004;f abbiano et al..20 04). 

The various observations fit naturally into a picture 
where outflows are caused by feedback from AGN accre- 
tion. This has motivated theoretical studies of the proper- 
ties of AGN winds (e.g., Shlosman, Vitello, & Shaviv 19851 
Naravan & Yi 1994; Konigl & Kaitie, 1994^.Stone & Normal 
19941 IMurravet alJ I1995t lElvisI I2000t IProga et alJ 1200(1 
Proga & jCallman f '200'4V The models have invoked vari- 
ous driving mechanisms near the black hole, including hy- 
dromagnetic disk winds (e.g., Konigl & Kartie 1994), MHD 
jet outflows (e.g., .Blandford & Pavne. .198 2) powered by 
magn etic coronae formation over the disk (jMiller & Stona 
120001), Comptonization in an X-ray halo (e.g., BegehnaM 
1985), and radiation pressure coupling to dust opacity (e.g., 
Dopita et al. 2002). Irrespective of their details, these pro- 
cesses universally yield a multi-temperature, clumpy, filled 
wind structure fBegelman & McKee 1983; Konigl & Kartj^ 
119941) which provides a good representation of the vari- 
ous c omponents and spectra l properties of observe d outflows 
(e.g., Krolik & Ki-iss 2001; Dopita et al. 2002; Ogle et alj 
^003). The winds evolve like a Sedov-Taylor type blast wave 
([Begelman et al. 198 3|), as suggested by observed velocity 
profiles (Shopbell & Bland-Hawthornll99 8). and develop the 
typical shell structure and phases of evolution of these blast 
waves (Schiano 1985h . well known from the analysis of su- 
pernova remnants (e.g.. [Ostr ik er & M cKee 1988). 

The theoretical works have examined the generation mech- 
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anisms of winds near a pre-existing accretion disk, and the 
observable impact of these winds on the in terstellar medium 
(ISM) as they shock and entrain gas (e.g.. 'Middelbe rg et alJ 
Il004; Machacek et al. 2004; O' Sullivan et al. 2005). How- 
ever, the AGN luminosity or accretion rate is generally an in- 
put parameter in these analyses, with an undetermined macro- 
scopic fueling mechanism. Without a self-consistent calcu- 
lation of the coupled evolution of the wind/blast wave and 
AGN accretion rate, such models, while critical for charac- 
terizing the detailed radiative properties of the ISM local to 
the AGN, cannot provide a physical motivation or understand- 
ing of the distribution of AGN luminosities, Eddington ratios, 
black hole masses, and fueling mechanisms, and their evolu- 
tion with redshift. 

If feedback from accretion couples to the gas surround- 
ing black holes in Seyferts, the observations and theoreti- 
cal models motivate the following picture for AGN activity 
in quiescent galaxies having a supply of cold gas. Obser- 
vations of the dynamics and distri bution of cold gas in the 
central regions of these AGN (e.g., Krolik & Begelman 
Kanekoetal. 1989; Heckmanetal. 1989; Israel et al. 1989- 
MeixneretaLlll 990; Granatoetal. 1997; Bocketal. 20od 
gghinnerer et al. 2000; Galliano et al. 2003; Radomski et aO 
2003; Weigelt et al. 2004; Jaffe et al. 2004; Prieto et al. 200i 
Elitzut 2005; M ason et al. 2005) suggest that rotationally sup- 
ported gas extends to the inner regions of the galaxy. Whether 
part of the galact ic disk or, potentially, a "clumpy" torus (e.g., 
lAntonuccilTl993t) or bar-like structure, molecular clouds or 
blobs of cold gas could be accreted stochastically by the black 
hole. Such gas will have some turbulent (random) motion and 
corresponding probability of colliding with the central black 
hole (see also e.g. lLauer e t al. 2005). 

The mass of gas required to sustain an AGN is only ~ 
10^ - 10^ Mo, far less than the ~ 10^ - lO'^M© observed 
within the inner tens to hundreds of pc of many late- 
type galaxies (e.g., Kaneko et al. 1989; Heckman et al. 19§^ 
iMeixner et al. 1990; Granato et al. 1997; Gal hano et alJ20ol 
iMason et aL .2005; .EUtzui .2005„ and references therein), 
a nd o nly ~ 1 event per Hubble time may be expected (see 
§ 14. 2> . Therefore, large-scale gravitational torques are not 
required, although mechanisms such as disk and bar insta- 
bihties (e.g ., Norman & Silk 1983; Norman & Scoville 1988; 
^hlosm an. Begelman. & Frank 1989; Lin. Pri n gle. & Reesi 
1^88; Lubow 1988), minor mergers (e.^.. [Ro osI 198i| 
Gaskell 1985; Hernauist 1989; Mihos & Hernquist' '1994 
^rnquist & Mihos 1995; De Robertis et al. 1998; Tanagucl^ 
19991) . or magnetic instabilities (e.g., Krolik & Meiksin 1990) 
may contribute to the "effective" turbulent motion. In a 
"collision," the black hole will accrete at a high rate for 
a brief period of time until feedback impacts the cold gas, 
driving a blast wave and initiating a feedback-dominated 
"blowout" phase. This blowout determines the subsequent, 
time-averaged evolution of the Seyfert light curve, obscured 
fractions, and Eddington ratio distributions, and the system 
decays to lower luminosities until a potential subsequent ex- 
citation. 

We use this picture to develop a model for the fueling of 
Seyfert galaxies, which allows us to predict their luminos- 
ity functions and evolution with redshift, among other quan- 
tities. In § 12] we derive a generalized Sedov-Taylor solu- 
tion for feedback-driven outflows and both a Bondi-Hoyle 
type and generalized perturbative accretion solution within 
such a medium. This is applicable to any feedback regu- 
lated accretion system, including the quiescent systems stud- 



ied here, "blowouts" in merger-driven quasars, and accreting 
black holes and neutron stars in supernova blast waves. In 
§|3lwe apply this solution to the cases of merger-driven activ- 
ity and stochastic accretion in quiescent systems. In §|3we 
calculate the expected excitation rates, duty cycles, and lumi- 
nosity function of Seyfert galaxies as triggered in quiescent 
galaxies. We further predict and compare with observations 
the distribution of host properties and morphological types as 
a function of luminosity, the distributions of accretion rates 
and Eddington ratios, and evolution (compared to e.g. evolu- 
tion in merger-driven quasar activity) with redshift. In §|5lwe 
predict the effects on the Mbh - c relation from this mode of 
AGN fueling, in particular estimating corrections to the slope 
and scatter at low Mbh- In §|6lwe determine the implica- 
tions of this picture for Seyfert obscuration and the potential 
buildup of the classical molecular torus. In §0we investi- 
gate effects on the host galaxy from this mode of fueling. In 
§111 we estimate and compare with other quiescent modes of 
AGN fueling, such as accretion of hot spheroid gas and stellar 
winds from passive evolution. Finally, in §|9lwe summarize 
and discuss predictions to test our model. 

2. ACCRETION IN FEEDBACK-DRIVEN OUTFLOWS 

A black hole accreting at the Eddington rate will grow ex- 
ponentially on a relatively short timescale (a few xlO^yr). 
If the black hole is small, feedback is suppressed exponen- 
tially, and surrounding gas can equilibrate with the low rate of 
energy or momentum input. However, if the black hole is al- 
ready large, it will be sufficiently luminous that the energy in- 
jected cannot be radiated by the gas in a local dynamical time, 
effectively resulting in an instantaneous, point-like (relative 
to the scales of the galaxy) energy injection in the center of 
an (approximately spherical) bulge which dominates the local 
gravitational potential. Therefore, it is appropriate to describe 
this phase as a Sedov-Taylor type blast wave, which, in detail, 
a number of simulations have shown is a surprisingly good ap- 
proximation to a full solution including radiative cooling, the 
pressure of the external medium, magnetic an d gravitational 
fields, and further effects, at least on galaxy jHopkins et alJ 
2006c) and larger ( Furlanetto & Loeb 2001) scales. 

Here, we determine the behavior of such a blast wave when 
energy is input from black hole accretion into a medium with 
an initial density gradient. In particular, we derive an ap- 
proximation to the internal density, velocity, and temperature 
structure of the blast wave as a function of time, which be- 
comes exact as r — > 0. Using this, we can calculate the evo- 
lution of the Bondi-Hoyle accretion rate for a central black 
hole at r = 0. However, because the Bondi rate is not calcu- 
lated self-consistently in a medium with external density and 
velocity gradients, we use our solution for the internal blast 
wave structure to derive a solution for a perturbative inflow 
driven by the black hole gravitational potential at small radii. 

Many authors have considered the evolution of blast waves 
and the detailed clumping and radiative properties of their 
interiors (e.g., [B egelman et al. 1983; Begelman & McKe^ 
'1983; Ostri ker & McK ee 1988, and references therein). (For 
a discussion o f the dif ferent Sedov-Taylor radial shells or 
phases see e.g. lSchianol lfr9 853.) However, these works have 
focused on understanding the radiative properties of such 
blast waves, and have not calculated an accretion rate solu- 
tion within them or examined the interior properties as they 
affect the small radii relevant for accretion processes. 

For simplicity, we assume a power-law scaling for the ex- 
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ternal (pre-shock) density, 

poWcx/^-*^ (1) 

appropriate for a spheroid or black hole-dominated potential 
or molecular cloud. Note that this need hold only over some 
range in r for each "stage" of the blast wave. We are, for now, 
interested in early times during the wave and (by definition of 
the blowout "triggering" condition) blast waves with energy 
greater than the binding energy of the material, so at least in 
this phase we can consider gravity to be a second-order effect 
which we calculate in detail below. 

Under these circumstances, the system will evolve as a sim- 
ilarity solution, in which the shock radius expands as 

7;,cxf''. (2) 
These similarity solutions are well-understood, and we refer 
to Ostriker & McKee ( 1988) for details. Once the system en- 
ters the blast wave-dominated phase, the central accretion rate 
will decline with time, in power-law fashion when the sur- 
rounding medium evolves according to a similarity solution. 
We derive this below; for now note that we expect 

L cx M cx f"''^ (3) 
with rji > 0. The total energy of the blast wave evolves as 



Es oc R 



-kE 



CX t 



VE 



(4) 



with rjE = -rjkE- Radiative losses from the blast wave 
front are small, at leas t on the scales of interest here 
iFurlan etto & Loeb 2001). However, if feedback energy con- 
tinues to couple to the surrounding medium as the accretion 
rate declines, there are two possible cases. First, if 

77i>l,7/£wO, (5) 
then the blast wave is dominated by the initial energy input, 
and the evolution is effectively energy conserving with kE = 0. 
Second, if 

0<riL<l, riE=l-TlL, (6) 

the energy input increases with time and radius in a self- 
similar manner In either case, we obtain a similarity solution 
for the shock radius of the form 

2 2 + 77£ 



77: 



5-kp + kE 



5 ka 



(7) 



2.1. Internal Blast Wave Structure at Small Radii 



We determine the internal structure of the blast wave using 
the two-power approximation (TPA) for the internal velocity 
and density structure. Unlike e.g. a one-power approximation, 
this ensures that gradients behave physically at small radii. In 
general, a quantity x will have the scaling 

x(r,t)=x,{R,)i{\) (8) 
internal to the blast wave, where 

\ = r/R, (9) 
and x\ is the post-shock v alue of x at the shock front. From 
lOstriker & McKeeUl988h . this gives 



p = apX^i''+{\-ap) A 



(10) 



where 



3-jk 



" = -(9-5^p) 



1 

7-1 2 
j(kp,cm-kp) 



5(2- kp) 



10-7-(7 + 2)fep 25+lUp 



Kp2 



:MIN(- 



-(^+l)(2^p-^ 



p, crit) 



epi,i)=n-8kp 



'^p.cnt ~ 



7+ 1 



= 2 



(11) 



(the second equalities are for 7 = 5/3). These are determined 
by requiring the gradients to match exact solutions at the cen- 
ter and shock front and imposing mass conservation. The cor- 
responding velocity approximation is 



, (7+l)A + (7-l)A^'^ 4 1, 
2^ -5^^5^ 



with 



42=1 + - 



7- 



p. cut 



-kp) = 6--kp. 



(12) 



(13) 



The expressions reduce to an exact solution for 7 = 5/3, 
kp = kp, ail = 2. For other cases, they provide a good approx- 
imation for our purposes, accurate to ~ 17% even for kp = 0. 
The accuracy improves at small radii, and these solutions cap- 
ture the correct power-law dependence as r ^ where, for ex- 
ample, Kahn's approximation ( Kahn 1976) for the density and 
veloc ity structure, which is accurate to ~ 4% JCox & Franco! 

reduces to the expressions above. These also provide 
a good description of obse rved blast wave velocity and tem - 
perature^tructure (e.g., Shoobell & Bland-Hawthorn 199« 
iCrenshaw et al. 200Q: Krplik & Ki'iss 2001; Smith & Wilsor) 
200U iKasniet alJ I2002t iVeilleux et alJ I2005^ iRiceetalJ 
2006tl . 

Given the density and velocity field, it is straightforward to 
determine the interior pressure. For an adiabatic blast wave, 
the entropy per unit mass is conserved in comoving coordi- 
nates, giving 

s = P/p^ = siM'^' 
M{r) = MM{X) 



"fkp-3-kE 
3-k, 



This yields 



f cx A 



VE 



(14) 



(15) 



for the pressure as f ^ 0. For a blast wave with no energy 
injection, then, we recover the well-known condition that the 
pressure is constant at the origin. 

In detail, we can use the equation of motion to solve for the 
pressure terms. Inserting the density and velocity fields into 
the equation of motion and solving yields a six-term power- 
law solution for the pressure 
p 

=■ =fl/>o + fl/>i X^'''^^ + ap2 X^"-'^^ 

+ flp3A^'"^^'^'+flP4A^''^^^-^' 

+ flf5A^'"^2^'+aP6A'''^^2'" (16) 



with 



api 



api ■■ 



api = 



apav,(l -a^ryi^i) 

(2 + ipi)r]Ui 
(ap-l)a^,(l-a,.?7t/i) 

(2 + £p2)riiyi 
ap(av-l)i-l+r]- 



A'77+flD(l+^l')??f^l) 



ap4 ■■ 



(l + ipi + e,)r]i^i 
(1 -flp)(a,.- + rj- £i,ri + av(l + £i)rivi) 



{l+lp2 + ty)r]vi 



ap5 = 



-flp(l -fl„)^4 



ap6 = 



(a, 



ipi+2e, 
-i)(i- 



lpi + 2l. 



(17) 
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where a,. = (7+ l)/(27) is the coefficient of the A term in v 
and I'l = v'l/vj = 2/(7+1). 

Because P = P\ at A = 1, we obtain a solution for opo- In 
terms of P 

(18) 



P = P(0) + flsA"'"+2 + 



this yields an expression for P(0) which is quite complicated, 
reducing for e.g. ^£ = 0, 7 = 5/3 to 



P(0)=^(82 + 



210 



266 



I'i-lQkp \9-Up 



7 



560 



108 



n-5ko U-5ko \6-lk. 



(19) 



which interpolates from P(Q) « 0.32 for /tp = to P(0) = 1 /25 
for kp = 2. However, the full expression for P(0) from the six- 
power expansion is reasonably well approximated (to ~ 10% 
for a variety of 7, kp, kE) by the simpler form 



P(0) = 



(7+l)'(^p. 



-kE) 



372 + 2O7+ 1 -(7+ 1)[(37+ l)^p- 3(7- 1)^£] 
AO-kp-kE) 



24-9kp + 6kE 



(20) 



which follows from the pressure-gradient approximation (see 
lOstriker & McKee 1988), essentially a two-power approxi- 
mation in the pressure. 

The next-order term, flpA^p'"*"^, which is proportional to the 
local pv^, is given by 



(7' - l)(fcp. cnt - kpm -kp)j-(2+7jE)) 

27(l + (2-fe,)7)(10-7-(2 + 7)^p)(2 + rte) 

8(2-fcp)(19-5^p-3??£) 
5(13-5fcp)(25-lUp)(2 + 77£)" 



(21) 



This goes to zero for j = 5/3, kp ^ 2, but this is because in 



this case a,- 



^p2 = 1 , and a proper solution yields 

2-riE 



2 + llE 



(22) 



It is straightforward to derive the remaining terms, but we are 
interested only in the behavior as r ^ 0, so we do not need 
higher-order terms in A. 

2.2. Evolution of the Bondi Approximation 

We can now determine the Bondi accretion rate as a func- 
tion of radius, again for r ^ 0, given by 



M 



(v2 + c2)3/2 



(23) 



where a is a constant of order unity dependent on the gas 
equation of state, p and = ^P / p ?at the local gas density 
and sound speed, respectively and v is the bulk motion of the 
gas relative to the black hole. Near the origin, P P(Q) = 
constant, thus 

PiO) = P(0)Pi=P(0)9oPoiRs)v^s 
7+1 

p^apX'^'"pi=apXiX'^'"po(Rs) 
7+ 1 



XI = 



7-1 



2pXi 



P 

— =Xi 



f-fP{m)\-^l^ ( r\¥.^ po(R,) 
\RsJ vARs)^ 



cx r 



^ flpXi 
epip(i-nyri(kp+iipi) 



(24) 



So, at some point, an inward-moving Bondi-like flow can 
grow and give a residual accretion rate, and if the shock wave 
continues to blow mass out of the central regions, the scaling 
should be similar to the time scaling above. The effective "ac- 
cretion radius" at which this will occur is, however, ambigu- 
ous. A natural choice is the black hole radius of influence, 

GMbh 



R 



BH ■ 



(25) 



where a is the bulge velocity dispersion. This follows from 
equating the black hole potential to the external potential, and 
determines where the black hole dominates the dynamics. We 
could instead consider the trans-sonic radius, 

_ GMbh . 

however, once Cs(/?bh) < cr, the trans-sonic radius will be 
larger than the black hole radius of influence, regardless of 
how C5 scales with R, and the black hole will not affect the 
dynamics. For the inner regions (or intermediate regions for 
kp = 2, 7 = 5/3), we have c, ~ v ~ r/t, so /?ts ^ ^bh at 
t Rbh/ a ^ ptdyn, where 

Mbh 

H = -p- « 0.001 (27) 

is the ratio of the black hole to the (total) spheroid stellar mass 
(e.g.,|Marconi & Hunt 2003) and 



a 7/(7 

fdyn EE - sa 3.6 X lO'yr -— — 

a \ lOOkms ' 



(28) 



is defined as the spheroid dynamical time (with a being the 
bulge scale length). Thus the timescale for Rts ^ Rbh is 
~ 10** yr, much less than the timescales of interest, and so it is 
inappropriate to consider accretion through 7?ts as it is outside 
the region where the black hole will have a significant effect 
on the kinematics. Even if some of the post-shock gas avoids 
cooling adiabatically and slows/compresses to pressure equi- 
librium, this will necessarily give ^ a, so /?ts ^ ^bh and 
again Rbh is the appropriate radius to consider 

Since we expect the accretion rate to decline rapidly with 
time, the blac k hol e mass is approximately constant over this 
period (see § 15.41 for a detailed calculation), so for a fixed 
accretion radius ~ Rbh this then gives 



p 

M oc — (X f 
ci 



3{l-V)-v(.kp+'jipO 



oc t 



(29) 



(30) 



This equation for r]i expands to 

_ 24-97-5fep _ 3(9-5;fcp) 

'^'■"(7-l)(5-^p)" 2(5-kp) 

for a rapidly declining accretion rate with r/£ = (i.e. valid for 
riE > I, kp < 17/13 for 7 = 5/3). For a less rapidly declining 
accretion rate with t/e y^O (i.e. < t^l < 1) we obtain a self- 
consistent solution with 

3(19-47-(4 + 7);tp) _ lll-5Up 



(16-5^p)7-l 



11 -25 ko 



(31) 



which is valid for all kp where 77^ < 1 (kp > 17/13 for 
7 = 5/3). Together, these form a continuous class of self- 
consistent solutions for all < ^p < 37/17. 
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2.3. Perturbative Calculation of the Accretion Rate 

Although the above derivation gives a reasonable solution 
for M and its sc aling with time, in agreement with simulations 
iHopkins et alJl2006ci,) . it is not strictly appropriate to adopt 
the Bondi-Hoyle rate at some radius, as the full blast wave 
solution represents an outflow. However, in the Sedov-Taylor 
solution, a given mass shell interior to the blast wave deceler- 
ates with time, and eventually the gravitational potential will 
dominate the dynamics. Shells with smaller initial r/Rs begin 
with smaller relative velocities and decelerate more quickly, 
and will successively "turn around" and fall back establish- 
ing a Bondi-Hoyle-like steady accretion flow. Therefore, we 
expect a mode of inflow in the inner regions which grows 
relative to the decaying blast wave solution, with the inflow 
smaller at early times when the blast wave is propagating 
rapidly (although the inflow may decay with time in an abso- 
lute sense). The gravitational term in the equation of motion is 
a perturbation, which will induce corresponding perturbations 
in the velocity and density. We can, therefore, describe this in- 
flow as a perturbation about the exact blast wave solution. We 
consider the density and velocity perturbations 6p ^ p and 
(5,, ^ v, and demand that the first order perturbation to the 
mass inflow rate 



5M = 5M(t) = 4-TTr^(SpV+p S,) 



(32) 



be constant in radius (although not necessarily in time). Our 
method and this assumption are identical to that used to de- 
rive the Bondi solution, but considered in a medium with a 
pre-existing large density and velocity field (and consequently 
first-order \M\ =47rr^p|v| > \dM\). 
The first-order continuity equation is 



86 



1 d 



-Ji+ SM(t) = 0, 



dt 



dr 



(33) 



so d6p/ dt = Q and Sp = 5p(r) because 5M is constant in r. The 
equation of motion 

dv\ dP 
dr / 



dv 



df 



(dv Ov\ 
■PKWt^'^r) 



dr ^ dr 



(34) 



(where (f>{r) is the gravitational potential and represents the 
driving first-order correction) is then completely determined, 
if we assume the perturbation is adiabatic; i.e. if 



5P Sp 



P 



apui) 
apayX\v 



(r->0). 



(35) 



The full equation of motion then reduces to the equation of 
motion for the perturbation 



6M 
A\n5M 



A,^t+Ar5pr^ =-pr^t^ — 



dlnf 



■ + r](ipi + kp)-i'o(ipi + l) 



\d\Ti6o 



d\nr \apa„xiiyi 



r](£pi+kp)+l 



^0 



+{-f-l)-ipl 2 

\apa„x\vf 



VQ = ayViri = -. 

1 



(36) 



2.3.1. Exact Isothermal Sphere or Wind Solution 



First, consider 7 = 5/3, rjE = 0, kp = k 



p, cnt ■ 



2, for which 



(37) 



the power-law approximations we have adopted and formulae 
above are an exact solution for the internal structure of the 
blast wave. Note that for this case (as indicated above) our a,, 
is inappropriate and the equations above should take 

2ri 1 
7+ 1 2 

In this case, Hp = £pi = I, rj = 2/3, and the time dependence 
of pr^t^d(j)/dr cancels completely, and since SM is a function 
only of time and 6p is a function only of r, it must be true that 
A,5Mt goes to zero or a constant. Then, both reduce to the 
same solution, namely 

dln^M . ^3^^ 



dlnf 

1 /d\n5p 
9 V dlnr 



-1 



+ 11 



^ _pt^d(p 
r dr 



/4po(ro)t(roy 



ro / 9r ' 

where ro is an arbitrary normalization radius. The equation 
for Sp has a general solution with Sp oc r"" , but this is not of 
interest as the central mass diverges strongly. This leaves the 
specific solution for Sp. Let us define 

dr a \aJ 

where e.g. k^=\ for the potential at short distances (domi- 
nated by the black hole) or large distances (when the galaxy 
potential is effectively a point-mass) and fe^ w for interme- 
diate distances in e.g. a Hernquist (1990) profile. This gives 



36 ~/a2f(ro> 



^(^)po(.o)(^) 
V rofl / \a) 

\ fHvn / \a 



-2 k J, 



ll-2k^ \ roa 

11-2^0 Vfdyn 



(40) 



where the second equality comes from setting the arbitrary 
ro = a and using f^yn = a/a. We show in §|3]that the charac- 
teristic blowout timescales f (a) are ^ fdyn as defined here, in- 
deed giving Sp-€ipal least for early times when the blast wave 
speed exceeds the escape velocity (such a condition guaran- 
tees t(a) < fdyn). 

This solution determines a self-consistent power-law slope 
for the decay of the accretion rate SM. The perturbation 
SM is independent of radius, whereas the blast wave outflow 
M (X r^pv (X r*t~^ falls at small r, meaning that the pertur- 
bation will eventually dominate at small r, as expected. The 
perturbations Sp and S^ both decline with radius, meaning that 
each is individually also important only at small radii, and 
vanish at large radii r ~ Rs where the blast wave structure is 
dominant. The decline of SM with time is less rapid than the 
decline of the blast wave M, reflecting the fact that as the blast 
wave expands, slows down, and cools, the gravitational poten- 
tial and accretion solution will increasingly dominate (albeit 
still with a declining accretion rate as gas continues to flow out 
of the central regions and be shocked to virial temperatures). 

This perturbative approach to calculating the accretion rate 
does not, however, determine the absolute normalization of 
that accretion rate, essentially an arbitrary parameter in the 
above equations, but within the radius where the steady ac- 
cretion flow dominates, a Bondi-like spherical flow should set 
in and detailed simulations suggest that the Bondi-estimated 
normalization is a reasonable a pproximation (.Hopkins et all 
l2006dlFurranetto & Loebl200l . 



Fueling Seyfert Galaxies 



7 



2.3.2. Extension to General Mass Profiles 

We now extend this solution to cases where kp ^kp^cm (i-C- 
7 ^5/3 or kp ^2). Because the power-law approximations to 
the density and velocity profiles are then not exact, the cancel- 
lation of the time dependence of pr^t'^dcj)/ dr is imperfect, but 
we can show that the remaining factor is small. If we neglect 
this small remainder, then a solution is obtained for A, = or 
AtSMt = constant. The latter case always yields the 



din SM 
dlnf 



■=-% = -! 



(41) 



solution. The former (A, = 0) yields one of two classes of solu- 
tions. If the accretion feedback decouples from the surround- 
ing medium after the initial injection (i.e. tje = 0, independent 
of tjl), then the general solution 



_ 2 _ 4 _ 12 

'^^"7''"7(5-^p)"5(5-^p) 



(42) 



is obtained. If the accretion feedback remains coupled to the 
blast wave, i.e. tje y^O for rji < 1, then a self-consistent solu- 
tion is obtained with 



18 



flL = 



2 + 7(5-/tp) 31-5kp 



(43) 



These are similar, with rji within 0.1 for a given kp, regard- 
less of whether the accretion feedback continues to couple to 
the surrounding blast wave. This implies that while such cou- 
pling may significantly affect the growth of the blast wave 
front, it does not dramatically change the time structure of the 
r ^ accretion rate once the initial blast wave has grown. 

For each, a corresponding solution for 5p is obtained, simi- 
lar to that for 7 = 5/3, jy^ = 0, fcp = 2. In all cases, the general 
solution for Sp has such a steep r dependence that it is gener- 
ally not of interest, but the specific solution is given by 



Sp cx r'''' 



^„i-l-2Av 



(44) 

which may increase with r, but always does so slower than the 
internal blast wave density profile oc r^^' and thus satisfies the 
necessary perturbative conditions. When A,6Mt =constant, 
i.e. ?7l = 1, a second specific solution also exists with 



Sp oc r 



(45) 



which may be of interest near the origin, and although the 
enclosed mass formally diverges, in reality the density profile 
will flatten once the perturbation dominates and a standard 
Bondi-type solution sets in, giving a Coulomb logarithm. 

We noted above that the time-dependent terms in 
pr^t^d(f)/dr do not entirely cancel in these approximations, 
and thus the solutions above are not exact. However, the 
residuals are small. If we apply the solutions above, then 
considering pr^t^d(f>/dr cx t^', we find that indeed S, is small, 
meaning this is not a bad approximation over any limited time 
range. The worst case is the tje y^Q, kp = I solution, for which 
S, = 3/8. However, generally the cancellation is more com- 
plete; for example for kp = Q we obtain (5, w 1/10, 1/20 for 
the rjE=0 and rjE^Q solutions, respectively. 

2.3.3. Late-Time Behavior of the Solution 

At late times, the internal blast wave velocity falls and the 
gas will not continue to cool indefinitely. The criterion for 
this can be determined exactly, but is essentially 



where t^y^ir) is the local dynamical time (distinct from the 
global dynamical time a/ a which we generally refer to). The 
entropy of the post-shock gas is large, and the sound speed is 
> (7, so as the gas cools and slows, a quasi-hydrostatic equi- 
librium is established, with 



dP 
dr 

„2 ^ 



1 dpc] ^ 
7 dr 
-(7- !)(/.. 



■PTr 



(47) 



If p can be separated into space and time dependent parts 
(as for a power-law dependence of p on time), the time de- 
pendence of p factors out of these equations, and remains 
nearly constant at fixed r while p{r) slowly declines as the 
residual outflow carries away mass. This approximation is 
valid (i.e. remains roughly constant given by the equations 
above while p decreases) so long as v ^ v'esc = \/-^M>- 

Then, we can approximate the residual velocity and decay 
of p with the original blast wave solution extrapolated to these 
late times. Strictly speaking, this is not appropriate, as at these 
times the potential becomes important and the solutions above 
are not necessarily valid. However, we can either: (1) con- 
sider special cases in which they remain a valid solution, or 
(2) perform a perturbative analysis again, but this time the 
blast wave solution is the perturbation to the quasi-static so- 
lution on large scales, determining a net decline in the density 
of the central regions. The perturbative analysis demonstrates 
that the blast wave solution is a self-consistent perturbation 
solution for r <C Rs, implying that the late-time behavior 



p 

M cx — r cx p cx f 

c 

is indeed appropriate. This gives 



'3-k 



(48) 



7- 



with ?7i > 1; i.e. rjE = for all self-consistent solutions, as 
should be true for the late stages of the accretion rate evolution 
when the blast wave is weak and expanding far from the ori- 
gin. Note that for the special case 7 = 5/3, fcp = A:p.ciit = 2, the 
proper solution becomes ryi = 2 (instead of 77^ = 1 as implied 
above, this is because the £pi is not valid in this specific case 
as noted above, and the exact £p = I should instead be used). 
Appropriately, this solution reduces to the scaling derived by 
IShu (.1977 ) for the self-similar collapse of a singular isother- 
mal sphere (kp = 2) with no initial velocity field, as the bulk 
outflow velocity is increasingly small relative to the escape 
velocity. This solution is important for late-time galactic- 
scale outflows, at low M, but is less relevant for e.g. accretion 
of cold clouds or supernova remnants, as the dense medium 
through which the blast wave propagates is sub-galactic. 

3. TIMESCALE FOR ACCRETION RATE DECAY 

3.1. Merger-Driven Black Hole Growth 

Now, we wish to determine the appropriate timescale for 
the accretion rate evolution calculated in § |2l It is convenient 
to employ the dimensionless accretion rate m = M/Mndd, with 

^Edd _ Mbu 
^2 



MEdd- 



(50) 



t > tdynir) 



(46) 



where e,. ^0.1 is the radiative efficiency and ts = 4.2 x 
10^(e,./0.1)yr is the Salpeter (1964) time. 

We first consider the case of gas inflows produced dur- 
ing a galaxy merger. The quasar light curves resulting from 
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this process have been discussed by Hopkins et al. (2005a-e; 
2006a-e); here we are specifically interested in the "blowout" 
phase. As the black hole mass increases exponentially in the 
late stages of a merger, a threshold is eventually crossed where 
feedback superheats the surrounding gas, unbinding the gas 
before it can cool. In ,Hopkins_et al ( 2Q06c ). we study this 
phase during galaxy mergers, and find that a similarity solu- 
tion with 

OTcxr''s (51) 

with a typical rji^l approximates well the accretion interior 
to the blast wave, ^opkinsetal. (200^ also find a weak 
dependence of -qi on the final black hole mass, expected from 
the differences in external profiles as derived in §|2] 

There are three regimes relevant to our calculation of black 
hole growth. First, there is a period of rapid (high-m) accre- 
tion, where accretion rate is high enough to be Eddington- 
limited. This proceeds efficiently as the black hole grows 
at low masses, radiating at low enough luminosity that the 
nearby gas can still cool and be accreted. 

The second regime is the beginning of the blowout - when 
the black hole becomes large enough that e.g. the feedback 
energy input into the surrounding medium in a local cooling 
or dynamical time is sufficient to unbind the gas. Effectively, 
the gas is "superheated" in the potential of the galaxy and 
the medium enters the Sedov-Taylor type solution of Equa- 
tion (in}- Of course, since the implied Bondi-Hoyle accretion 
rate is much larger than the Eddington rate, it will take some 
time tm for the Bondi rate to fall to MBondi < A/edd- 

Beyond this point, the third regime is entered and the actual 
accretion rate begins to decline. If f = at the beginning of 
the "blowout" (second phase) and f„-, is defined as the time for 
the Bondi rate to fall to mBondi = m= 1, then the subsequent 
accretion rate is given by 



t \-nL 
m = \ — 



(52) 



We did not attempt t o calculate the timescale t,,-, from our ana- 
lytical model in Hopkins et al.i (20Q^, as it does not affect 
our analysis of the faint-end slope of the quasar luminos- 
ity function. However, this quantity determines the relative 
contributions of e.g. stochastic accretion and relaxation from 
bright, high-;7i phases in mergers, so we calculate it here. 

During the blowout, the accretion rate drops to that given 
by the internal structure of the blast wave in Equation (I24> 



(y2 + c?)3/2 



(53) 



= 47rax 



-3/2 



The relevant radius for calculating the accretion rate is, as 
discussed in § |2 the black hole radius of influence Rsh = 
GMbhO'~^- For convenience, we normalize the accretion rate 
to Rs = Rbh and define Jbh = t(Rs = Rbh)- 
The accretion rate for f > fee is then m = mo (f /tsuT^'', with 



/7P(O)0o\-3/2 G^MBHPots 



(54) 



apXi ' v,(RB\\y' 
and r\i^ is as in Equations (I42> and (I43> . Therefore m = 1 at 

fm = fBHmi^'''. (55) 



To solve for nio, we consider the properties of the central 
regions of the merger, for which the stellar bulge has, in gen- 
eral, nearly finished forming by the time the "b lowout" oc- 
curs (Sprinael et al. 2005a; Hopkins et aLll20063) . We define 
a characteristic bulge radius a by 



GMbu 



(56) 



which is similar to the half-mass or half-light bulge radius 
(with the exact proportionality constant relating the two de- 
pending on the shape of the bulge profile). The definition 
for 7?BH gives the relation 7?bh = p.a (where /i is the ratio of 
black hole to bulge stellar mass). We can relate the spheroid 
and blac k hole properties v ia the black hole-bulge mass re- 
lation of Marconi & Hunt (2003) and Mbh - o" relation of 
iTremaine et aL C2002) (see Equation 1891 ). which then imply 



a = flo I — 
Vcro 



yn ■ 



a 



(57) 



(58) 



where we have defined the bulge dynamical time above and 
ao = 14.5kpc, fdyn.o = 7.13 x 10'' yr for ctq = 200kms"'. It is 
also useful to define an effective gas mass fraction (just before 
the "blowout") 



Po(^Bh) = /gas 



3Mbui [Rbh 



(59) 



47ra-^ \ a 

where /gas is the total bulge mass fraction in gas (/g^s = 
Mgas/Mbui), evaluated at the time of the "blowout" (usually 
^ 1 %) and p* is a correction for the shape of the galaxy profile 
(i.e. for the fact that our definition of a is not equivalent to the 
half mass radius), with e.g. p* = [2(2 + V 2)2l/[3(3 + V2)^1 = 
0.07 and p* = 1/(3 tt) = 0. 1 1 for lHernauisK 1990.) and isother- 
mal sphere profiles, respectively. This gives 



G/Oo(^bh) = 



/gasP* 3 fR 



'dyn 



Combining these equations and scaling by cja gives 
ts_\ /c,,o(^bh)\"V^bh^"*° 

^ f dyn ' 



mo = 3a/gaspV(- 



)(^BH) y^^^BHy 



(60) 



(61) 



The energy condition for blowout, namely that sufficient 
energy be released to rapidly unbind the surrounding gas, is 

given by Vs(/?bh) = I'esc, with a well-defined Vesc = V2cr at 7?bh- 
Therefore, with these choices, 

fBH = V—TZ T = ^M^dyn ■ (62) 

We can then use Equation ( I55t to determine 1^: 

^. = .(3a/ga.,P*)'/-(|^)'"''"(^ 

^i.(i-..)M(£M^^)-/'^^, (63) 

For the simplest, self-similar example of Sedov-Taylor ex- 
pansion under conditions resembling our simulations, we 
have rji = 2, rj = 2/3, which corresponds to the solution in 
e.g. an isothermal sphere, wind, or any v cx r velocity struc- 
ture, and indeed this provides a good first-order approxima- 
tion to the blowout p hase over a wide range of conditions 
jHopkins et alJl2006cl) . In this case, we can simplify the 



Fueling Seyfert Galaxies 



9 



above equation. We expect Vs(/?bh) ^ o", as is essentially 
guaranteed by the blowout condition that the sudden energy 
input be sufficient to unbind the gas locally. We also expect 
c.5.o(^bh) ^ 10 kms"', since it is inflows of cold gas from the 
merging disks which are being accreted. This gives 



f„i=4.08 X lO^yr f 

/ /gas 



Mr 



1/2 



.0.01 



1/2 



(64) 



0.001 



Q,(,(-Rbh)\-3/2 



10km s- 



Here, we have again used the M^h - c relation to replace the 
<T dependence with an Mbh term. Direct comparison with 
the best-fit to the blowout in our simulations shows such a 
trend, with a scatter of a factor ^ 2 about the relation above. 
This scatter is smaller than that implied by the various factors 
above, but this may be because several are correlated. 

3.2. Molecular Cloud Accretion 

Next, consider the accretion of a molecular cloud. As we 
describe in § 14.11 we expect clouds to be moving at a speed 
~ lOkms"' near the black hole, with characteristic gas densi- 
ties n ^ lOOcm"^. This gives an initial Bondi-Hoyle accretion 
rate m ^ 1 (limited to m = 1) for all black hole masses of in- 
terest. Provided that the black hole has grown large enough 
to lie on the Mbh - cr relation (see § |5}, it will release suffi- 
cient feedback to unbind the molecular cloud in a timescale 
< 10^ yr. Because this is small compared to the Salpeter time, 
the cloud crossing time, the cloud dynamical time, and the 
timescale for the moderate and low-accretion rate phases of 
the "blowout", we can effectively treat such an event as if 
it enters the blowout phase immediately upon the "collision" 
with the cloud. 

The calculation of the relevant timescale for the accretion 
rate decay is identical to that in § 13.11 since the typical cloud 
radius is larger than the black hole radius of influence, espe- 
cially at the relatively small masses of interest. However, the 
surrounding medium is not the same, and the condition for 
"blowout" is modified by the need to unbind a local cloud of 
cold gas as opposed to the entire gas content of central re- 
gions of the galaxy. For a molecular cloud, we expect kp = Q 
(at least effectively, since the black hole does not necessarily 
lie at the cloud center), -q = 2/5, and po{Rbh) = Po = 'WpWci 
with Hci ^ lOOcm"^. The speed at some radius and time 
is given by v, = rjR^ii/t. We again define Jbh as the time at 
which Rs = ^BH- 

The similarity condition for the blast wave gives 



^Eb(Rs) 



p{Rs) 



1/5 



.2/5 



(65) 



where Ef, is the blast wave energy, p = po is the mean density 
inside the blast wave (necessarily equal to the mean density 
outside by mass conservation and since kp = 0), and 

^=ir^T- (66) 



477-1]- (7 h 

= Eh/Mv] 



are generally constant for a self-similar blast wave (and ab~l 
for the early stages of the blast wave where Vj ^ Vesc)- 

Consider first the simple case of an energy-conserving blast 
wave, in which the accretion rate declines quickly enough that 
the energy input is dominated by the initial accretion (or ap- 
pUcable if e.g. the ionization or initial blowout of the inner 



regions de-couples the blast front from the black hole feed- 
back). Simulations of particular wi nd mechanisms find tha t 
driving is efficient for m > 0.1 (e.g.. lBalsara & Krolikll993h . 
which also suggests this is an accurate approximation. The 
blowout condition essentially defines the blast wave energy, 
as that required to unbind and expel the cloud 



Eb ^ -^binding ~ Md (T^. 

This implies rj = 2/5, and plugging this in we obtain 

R. ,Xl2(Rcl\^l^ 



where 



Rs A 12 ( RciY 



fb^(^^ 

^ ^bindina 



(67) 



(68) 



(69) 



is approximately constant, at least until the late stages of the 
blast wave evolution. 

Using this result in Equation \5A\ . again normalizing to 
7?BH, we obtain an equation for rii{t). We again re-normalize 
to m = {t/t,i,)~'^<-, giving t,f, = Jbh '«(>'''"■ ■ Solving gives 



_ ^BH Tj A/jjifRsH 



3(r,/'-l)/2 



mo ; 



47rQXi / lP(0)eo \ g2 Mbh Pq ts 



(70) 



Finally, we use kp = 0, tje = 0, adopt the Mbh - c relation as 
an initial condition, and take 7 = 5/3 (a = 1/4 for the Bondi 
solution) to obtain 



f„i = 0.40 — n 
cr 



(71) 



mo 



^•^^^"'Hloo^ 

cr 

lOOkms- 



lOOkms-i. 

* VlOOkms-iy VlOOcm-3/ VlOOpc/ 



where the t* term is a small, order unity correction which 
depends weakly on the properties of the system. 

This more detailed analysis recovers the expected scaling, 
that the characteristic timescale for accretion decay is given 
by the cloud size divided by the local escape velocity; i.e. 
approximately the velocity at which the cloud will be unbound 
and expelled. Although we have solved this in detail for the 
7 = 5/3, A;p = case, a similar calculation for different values 
of 7, kp, rji, and rjE gives a nearly identical equation for t,n, 
with small, order unity numerical changes to the coefficient 
of /?ci /cr as expected from Equation MQ\ and slightly different 
power-law coefficients for the terms in t*, with similar weak 
dependencies on cr, no, and R^i in . 

One important point to check is that our approximation of a 
constant Mbh is reasonable. For f„i we have 



t,i, w 3.9 r, X 10^ yr ( ] , 

^ VlOOpc/ VlOOkms-1 



(72) 



which is much less than the Salpeter time (f,,-, ^ 10"^ ts), and so 
the black hole mass is roughly constant during blowout. This 
also means that once the system enters blowout, the accretion 
rate drops below m = 1 immediately, and so a pure-power law 
decay is a good approximation when considering the proba- 
bility of viewing the black hole at a given m. 
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For the luminosity-dependent Seyfert lifetime, a power-law 
accretion time history of the form in Equation (I52> gives 



df 



dlogm 



: Tjitii, lnlO;7i' 



-UVL 



(73) 



which we can use to infer the probability of observing a sys- 
tem at a given accretion rate. Since Mbh is roughly constant 
once the blowout begins, this can be directly translated to a 
distribution in luminosity using L = riiL-EM = merMButs^ c^- 

4. ACCRETION OF MOLECULAR CLOUDS 

4.L Rate of Cloud-Collision Events 

If we consider a black hole in a system with a significant 
disk or other cold gaseous component which extends near the 
center, then there will be some probability for the black hole 
to cross paths with a dense molecular cloud. The timescale 
between collisions with clouds is 



revent = (n7r/?^lVcl) ' 



(74) 



where n = Nd/V is the number density of cold clouds, R = 
Rci is their typical size (giving a cloud cross-section tt/?^]), 
and v'ci is their characteristic random velocity. Although the 
gravitational interaction with the black hole is a long-range 
force, the effective cross section for interaction with a cloud 
will be ; 



7r/?^[ if the black hole radius of influence. 



^BH — 



GMu 



Mr 



lO^M,, 



1/2 



(75) 



is much smaller than the cloud size, which for a typical large 
molecular cloud with R^ ^ 10- lOOpc is satisfied for all black 
hole masses of interest, especially for the small Mbh charac- 
teristic of late-type systems. 

The typical random cloud velocity is essentially the vertical 
velocity dispersion of the disk, Vd = cf'^^ ^ lOkms"^ in order 
to maintain pressure support. In any case, the value of cf'^^ 
is relatively unimportant as we demonstrate that it ultimately 
cancels below. The cloud number density is 



"ff 



N 1 A^Vci 



V Vol V (4Tr/3)R 



(76) 



where v// is the volume filling factor of cold clouds, so 



1 Rc 



.disk 



10"' yr 



VO.OOl/ VIOOdc/V 



„disk 



lOOpc/ VlOkms- 



(77) 

Therefore, a large fraction of black holes in disk-dominated 
systems should have undergone such an event, although we 
determine the effective "duty cycle" more completely below. 

Next, we consider how vyy and fevent scale with host galaxy 
properties. Because most of the gas in the ISM resides in cold 
clouds, the cloud filling factor is 



PiSM "ISM 



Pel 



"cl 



(78) 



(where pisM and pc\ are the average ISM and cloud gas den- 
sities, respectively) independent of the size distribution or 
shapes of clouds. For an exponential disk with total gas mass 
Mgas = /gasMd, the surface density profile is 



-gas('-) = -:; — 2~ exp(-r/rd) . 



(79) 



In the central regions of the disk (taking r 0, z 0), this 
resembles an isothermal sheet with constant scale height h, 



PiSM ■ 



47rr^/i ' 



(80) 



The scale height h is in general given hy h ^ rd{cf^^/vc), 
where vv is the circular velocity. As r ^ 0, the potential 
is dominated by the bulge, i.e. Vc(disk) while the bulge 
a constant, so h ^ (cf^^/ a), giving 



PlSM 



a 

^disk ' 



(81) 



For a disk formed by collapse at fixed baryon fraction and spin 
parameter A, Md oc and the quantity 



1 Md 1 Md 



lip Anrl 



i 0.4cm 



MW 



(82) 



is constant across disks, where we take M,/ = 10" M© and 
Rd = 9kpc for the Milky Way in normalizing this; i.e. essen- 
tially a baryon fraction « 0.05 and spin parameter A = 0.04 
for Vvir w 160kms"'. Note that for this Vvir, calculating the 
gas density with /gas = 0.1 at r = 8kpc (including the scaling 
of circular velocity, scale height, and surface density to this 
radius) gives the standard «ism ~ 0.3 cm"^ for the local ISM. 

For a standard cloud density Wci = lOOcm""*, this then gives 
for the volume filling fraction at the center of the disk 



where 



pel 



Md 



(83) 



nipn^i 4-nrl 



i4.0x 10" 



«cl 



100 cm- 



(84) 



is approximately constant. The characteristic timescale for a 
collision with a molecular cloud is then 



1 Rc 



/gas" cr 



(85) 



and the factors of cf^^ cancel (as a higher cf^^ will increase 
the speed of clouds and thus number of collisions, but also 
"puff up" the disk and reduce the density of cold clouds). In 
addition to producing an appropriate Eddington ratio distribu- 
tion, there is observational support for stochastic cloud colli- 
sions with AGN on such a timescale, as e.g. dust clouds in the 
central regions of nearly low-luminosity AGN appear to trace 
accretion events with a similar frequency (Lauer et al. 2005). 

The expected spectrum of cloud sizes is n{R R + dR) cx 
R~* dR so that the filling factor for clouds in some interval of 
R goes as dR/R (xdlogR oc c/logMd, and there are equal con- 
tributions to V// from each logarithmic interval in cloud size 
from the smallest ^ Ipc to the largest > lOOpc. Therefore, a 
more complete calculation of fevent considering the probabil- 
ity of collision with clouds in each size interval gives only a 
logarithmic ln(/?ci,max/^ci,inin) correction to the rate at which 
clouds collide with the black hole and ultimate the "duty cy- 
cle" of activity. We subsume this factor into n, as it is within 
the present uncertainties. However, this range of cloud masses 
can be important for e.g. the scatter in the expected Mbh - c 
relation from events fueled in this fashion. 
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Finally, the mass of a clo ud needed for an accretion "event" 
is quite modest (see § l5.4t . as low as ~ lO^M©, and the ul- 
timate duty cycle derived in § 14.21 is effectively independent 
of the mass of cold gas inflow. Therefore, only a small frac- 
tion of the ~ 10^- 10^ M0 of cold molecular gas typically 
concentrate d within the central regions of late-type galax- 
ies (Kanek o et alJll989t iHeckman et alJll989 : Meixn er et al. 
[1990; Granato et al.'l 19971 iGalliano et alJl2003^ .Mason et al. 
|2P05; Elitzur 2005) needs to pass randomly near the black 
hole, and from our derivation above events are expected only 
of order once per Hubble time. With these considerations, 
there is no "angular momentum problem," at least on the 
scales of interest. 

Unlike bright quasars, which require large gravita- 
tional torques to sustain high accretion rates m ^ 1 
for - lO^-lO'^yr with large Mbh > 10** M© masses, 
no disturbance to the galactic gas, even in the cen- 
tral disk, is required for these quiescent Seyferts. Of 
course, pr ocesses such as such as disk and bar instabil- 
ities (e.g., Norman & Silk 1983; Norman & Scoville '1983 
Shlosman, Begelman, & Frank 1989; Lin, Pringle, & Reep 
1988; Lubow 1988), minor m ergers (e.g., |R oos 1981; 
Gaskell. .19 8.5: Hernauist 198"^, "Mihos & Hernquist 1994; 
Hernauist & Mihos 1995; De Robertis et al. 1998; Tanaguchi 
or magnetic instabilities (e.g.. iKrolik & Meiksin 199Q() 
may nevertheless play a role, and it is still of interest to un- 
derstand in detail the transport of material from scales of order 
tens of pc considered here to the small ^AU sizes of an accre- 
tion disk, but these processes, at least on the large scales we 
consider here, will ultimately serve to increase the effective 
random (non-rotational) velocity dispersion of clouds. 

There are however two limits in which this activity will be 
suppressed, i.e. fevent can be become much larger than a Hub- 
ble time. The first, in which /gas ^ and there is simply no 
gas supply to fuel this mode of accretion, is straightforward. 
The second, in which cr ^ 0, implies that the central gas in the 
AGN is not dynamically hot - i.e. it has relatively little dis- 
ordered motion and there is no disturbance which can bring 
gas to the black hole. Such systems will also have very small 
black holes, given the Mbh - c relation, so their contribution 
to t he ob served Seyfert luminosity function (which we discuss 
in j4!4l will be small. 



4.2. Light Curves and Duty Cycles 

When colliding with a molecular cloud, a black hole on the 
A^BH - relation will immediately enter the "blowout" phase. 
Even if the black hole is slightly undermassive and accretes 
before entering this phase, this time is short compared to the 
eventual tim e at low or m oderate accretion rates (see the dis- 
cussion in § l5.1l and l5.3t . and the distribution of duty cycles 
can be well approximated by neglecting these times. In § 13.21 
we determined that the light curve in such an event is given 
by m = (r ItiiiY"^^ (Equation 1521 ). where f,,-, cx Rc\/(j is given in 
Equation Mil . This yields a differential time per logarithmic 
interval in luminosity, dt/d\ogm oc t„;m~^/'^'', given in Equa- 
tion il3\ . The effective "duty cycle" for an object as a func- 
tion of the accretion rate m is then 

d/ _ dWevent d/ 

dlogm df dlogm 

= r]L— In 10 m-'/"' = So m"'/"^ . (86) 

^event 



Using our solutions for fevent (§ 14. U and t,,-, (§ 13. 2t . we obtain 
So ~ 10-' (^^)VLT, InlO, (87) 

where is defined in Equation i72\ . 

In general, 77^ ~ 1, so we can easily infer some important 
properties of the duty cycle. First, this implies a duty cycle 
at large accretion rates m > 0.1 of 1%, similar to that es- 
timated observationally from e^g. Kauffmann et al. (2003a); 
lYu ef a l. (2005); Do ng & De Rober tis (2005). The duty cy- 
cle becomes large at lower accretion rates, going to ~ 1 at 
m < 10"^, again simi l ar to that measured observationally from 
e.g. lHao et all ( l2005l) : lBest et alJ ( 120051) who find a large frac- 
tion of late-type galaxies hosting moderate/low accretion rate 
Seyferts. Note, however, that technically these are theoretical 
upper limits to the duty cycles, for if accretion proceeds inter- 
mittently (i.e. in short, potentially super-Eddington "bursts") 
the same average accretion rate on the timescales relevant for 
our calculations is maintained, although the timescale for such 
bursts is still constrained by the observed episodic quasar life- 
time (see e.g., M artinL200 4). 

This also implies an effective minimum accretion rate Jiimin, 
at which the total time spent with m > nin^in is equal to fevent- 
This is not to say that this is a hard minimum to the Seyfert ac- 
cretion rate, but rather that by the long timescales expected for 
decay to accretion rates much below this, a subsequent colli- 
sion with a molecular cloud is expected, which will "reset" the 
system to a high accretion rate. In detail, for m < m„im with 
Poisson statistics for the excitation rates, an exponential cut- 
off is expected to introduce a term exp [-(df /dlogm) /fevent] ■ 
Because this is a rapid cutoff, we can temporarily approxi- 
mate it as an absolute cutoff and determine rrirmn as where the 
duty cycle / = J d/ 1 . This gives 

which for < 1 gives mmin ^ 10""^. 

Finally, note that the terms involving cloud sizes have com- 
pletely canceled in this derivation. Thus, while the cloud size 
may be important for e.g. scatter in the Mbh-c relation (see 
§ 13, it does not enter into our ultimate calculation of the 
Seyfert luminosity function and distribution of accretion rates. 
Therefore, the considerable uncertainties in the properties of 
clouds at the centers of galaxies, and the "typical" giant cloud 
size related to Seyfert activity (and the immediate source of 
such clouds) do not affect our calculations. Rather it is the 
well-constrained quantities such as gas fraction and our theo- 
retically determined rji which determine these predictions. 

4.3. Seyfert Luminosity Function 

From the duty cycle as a function of accretion rate, we 
can determine the Seyfert luminosity function implied by 
this mode of fueling. We assume the black hole mass 
remains roughly constant during the blowout (see § \5.4l . 
giving d//dlogL = d//dlogm with L = merM^n/tsc^ , and 
assume a constant Cr = 0.1, as expected for a standard 
.Shakura & SunvaevI ( fl97l) thin disk. 

To derive the expected Seyfert luminosity function, then, 
we require the distribution of black hole masses and corre- 
sponding host galaxy properties, in particular their masses, 
gas fractions, and velocity dispersions. We use the lo- 
cal galaxy lu minosity functions from the CfA redshift sur- 
vey in B-band dMarzke et alJll994altJ) separately determined 
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for each morphological classification of E, SO, Sa/b, Sc/d, 
and Sm/Im. We convert these to ^"20 luminosity func- 
tions with the conversions for each type from this sur- 
vey as in Kochanek et al. (2001) (see their Table 5), and 
then to mass functions with the typical TT-ban d mass-to- 
light estimates based on mass and type from iBell et all 
( 120031) . which incorporate corrections for galaxy evolution 
and detailed stellar population synthesis from the models of 
IFioc & Rocca-Volmerange ( 1997). 

Grouping the E and SO galaxies as "early-type" and Sa/b, 
Sc/d, and Sm/lm galaxies as "late-type," we compare these di- 
rectly to the K-hand luminosity functions of Kochanek et aTl 
(l2001h . and find agreement (as the authors derive). We sim- 
ilarly compare to the m ass functions of early and late-type 
galaxies determined by Bell et all ( 1200 3) in both g and K 
ban d, and find reasonable agreement. Because the more re- 
cent |BelOIH] |2003) mass functions are estimated from the 
much larger combined SDSS-2MASS local galaxy sample, 
with more detailed correction for stellar mass-to-light ratios, 
we re-normalize our mass functions slightly (a small < 30% 
correction) to reproduce the Bell et al. (2003) late and early 
type mass functions. 

Next, we estimate the bulge-to-disk ratio of each morpho- 
logical type following Aller & Richstone (2002); Hunt et aQ 
( 12004a ) in the Z?-band and //-band, respectively (adopt- 
ing a similar procedure to correct to a mass ratio). This 
gives a bulge-to-total mass ratio of approximately B/T = 
(1.0, 0.57, 0.19, 0.077) [±(-, 0.2, 0.4, 0.5) dex] for galaxies 
of types (E, SO, Sa/b, Sb/c), respectively. The bulge-to-total 
mass ratio of Sm/lm galaxies is uncertain, and these galaxies 
may have no bulges whatsoever, but in any case we find be- 
low that their contribution is sufficiently small that they can 
be neglected even taking a maximal B/T for such systems. 

From the stellar mass function and B/T ratio for each mor- 
phological category, we construct a cumulative bulge or disk 
mass function and mean B/T as a function of mass. We 
do so and compar e with the bulge and disk mass functions 
of lTasca &Whit e' (12005), and the estimated mean B/T as a 
functio n of mas s from the size and surface brightness anal- 
ysis in IShen et aL (i2003i) . and find agreement in both cases, 
suggesting that this decomposition is reasonable. Although 
more detailed properties such as the mean disk gas fraction 
/gas are generally unnecessary for our subs equent analysis, we 
determine them from the compilatio ns of iRoberts & Havnea 
( fT99l andlKauffmann et al.' ("200313!). 

From th e bulge mass function, t he MR H-M hnire relation 
jMaaorrian et al. 1998; Marconi & Hunt 2003: Hiring & Rixl 
12004) then determines the black hole mass function. We 
adopt the relation of Marconi & Hunt (2003), which is es- 
sentially identical to that measured in simulations of spheroid 
and black hole formation (Robertson et al. 2005), and is also 
equivalent to the observed Mbh - c relation determined in 
iTremaineTt al. (.2002) , gi ven the relation be tween radius and 
stellar mass determined in lShen et alJ (l2003h . This gives 

Mbh = 0.001 Mbuige 

^- = ^«'"^«(200kJi^) • ^''^ 
The scatter in these relations is observed (and determined 
in simulations) to be » 0.35 -0.4 and 0.3, dex, respectively 
(see also Novak et al] ( 12005ft ). We assume that the PDF for 
black hole mass at a given bulge mass is distributed lognor- 
mally, with a dispersion equal to the above dispersions. For 
each bulge mass, then, we convolve the bulge mass function 
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Fig. 1. — Inferred black hole mass function from the moiphological 
mass f unctions of Marzke et al. 1 1994a b); Kochanek et al. i2()01); Bell et aO 
tlOOJ), with the typical bulk-disk decompositions of Aller & Richstong 
(2002); Shenetal. |2003); Hunt et al. (2004a); Tasca & White |2005) and 
MflH-Mbulse relation of Marconi & Hunt (2003); Robertson et aU >2005i) or, 
equivalently, Mbh — relation from Tremaine et al. 1 2002). The mass func- 
tion from E (red), SO (yellow), Sa/b (green), Sc/d (cyan), Sm/lm (blue; as- 
sumes a maximal bulge B/T= 0.02), and all types (black solid) are shown, and 
compared to the la range of the inferred mass function from Marconi et aU 
1 2004), based on the obs ervations of Marzke et al. 1 1994a); Kochanek et aU 
12001); Nakamura et aU ^2003); Bernardi et ah, i .20O3) ; Sheth et al. i200a 
(yellow shaded range). The results adopting the modified Mbh — O" relation 
suggested in §|5]are also shown (black dashed). 



with the PDF for black hole mass, and determine the resulting 
black hole mass function (BHMF). 

In Figure ^ we show the resulting BHMF, calculated 
individually for each morphological category of E (red), 
SO (yellow), Sa/b (green), Sc/d (cyan), and Sm/lm (blue). 
The cumulative BHMF obtained by summing these con- 
tributions is shown as th e black line. For c omparison, 
the BHMF determined in iMarconi et in (l2()0l is shown 
as the shaded yellow range, which a lso agrees wi t h othe r 
measurements of the BHMF by e.g. Salucci et al. (19991): 
Marconi & Salvati (2002); Yu & Tremaine (2002); Ferrare^ 
(2002); Aller & Richstone (2002); Shankaret al. (2004). That 
the agreement with our estimate is good suggests that our 
adopted conversions and decompositions are reasonable. For 
the Sm/lm case, we have assumed a maximal B/T = 0.02, 
only a factor of w 3 below that of the Sc/d galaxies, and it 
is apparent that the contribution to the integrated BHMF and 
number density at any mass of interest is negligible. There- 
fore, we ignore these galaxies in our subsequent analysis, as 
their B/r is in detail quite uncertain (if it is non-zero at all). 

We now have the observed BHMF, with the corresponding 
gas fraction, bulge size and velocity dispersion, and disk mass 
and size for each system. Our calculations above then allow 
us to estimate the rates, duty cycles, and lifetimes of Seyfert 
activity fueled by the quiescent accretion of cold gas, and the 
corresponding Seyfert luminosity function. 

For each black hole (i.e. each point in the joint PDF of 
■^BH, cr, /gas, Md, and host galaxy type), we convolve the 
black hole distribution with the duty cycle as a function of 
these properties. In other words, the observed luminosity 
function is given by 

0(L) = = I -^{x)dn{x) (90) 

dlogL J dlogL 

where x = (Mbh, /gas, o", Md), for example. Again, we have 
L oc m so dlogL = d\ogm. Given our blast wave solution 
for d//dlogm, the distribution of host galaxy properties then 



Fueling Seyfert Galaxies 



13 



completely determines the predicted Seyfert luminosity func- 
tion (insofar as it is attributable to this fueling mechanism). 

In Figure 12] we show the Seyfert luminosity function pre- 
dicted by our blast wave model and the assumption of fu- 
eling by the accretion of molecular clouds. Our predic- 
tion, in terms of the bolometric luminosity, is converted to 
a B-band luminosity function u sing bolometric corrections 
following [Marconi et al.' ^2004"), based on optical through 
hard X-ray observations (e.g., Elvis et al. 1994; George et al. 
[1998; Vanden Berk et al. 2001; Perola et al. 2002; Telf er et alJ 
12002; Ueda et al. 2003; VignaHetal. 2003), with an X- 
ray reflection co mponent generated by the PEXRAV model 
iMagdziarz & Zdziarskiiil995.) . and the appropriate Jacobian 
factors inserted. We have not calculated the effects of obscu- 
ration, and thus the luminosity function plotted is in terms of 
the intrinsic B-band luminosity, not necessarily that observed. 

For comparison, we plot the local (0 < z < 0.15 ) AGN lu- 
minosity functions from the SDSS, determined in lHao et alJ 
(|2005). The luminosity function estimated from Ha 
emission-line measurements is shown as blue points, over 
the luminosity range for which measurements exist. We 
show both the best-fit Schechter function (circles) and dou- 
ble power-law (squares) to the data. We convert to the B- 
band following Hao et al. (2005), based on the color correc- 
tions determined therein and from Schneider et al.'(e.g.,'200J, 
12003) . To show the effects of different bolometric corrections, 
the open square s adopt th e constant bolometric corrections of 
iMulchaev et alJ ( 11994.) : .Elvis et al.kl99 4). In either case, the 
authors note that this gives agreement with previous, but much 
shallower and more poorly constrained B-band AGN luminos- 
ity function measurements from Huchra & Burg ( 1992) and 
[Ulvestad & Ho (2001). The luminosity function shown is de- 
termined from the Ha narrow-line component, and is thus ex- 
pected to be isotropic to Seyfert 1 and Seyfert 2 galaxies, i.e. 
tracing the intrinsic luminosity which our prediction shows. 

We also show (red points) the luminosity function de- 
termined from [Oil] emission-line measurements, which is 
also expected to trace the intrinsic luminosity, with circles 
and squares again showing the best-fit Schechter and double 
power-law functions over the observed range. The correc- 
tion to the B-band foflows Mulchaev et al. ( 1994); Elvis et d] 
([1994), and gives agreement with the Ha determination over 
the range where the observations are most well-constrained. 
However, this does illustrate the potential importance of sys- 
tematic effects, especially at low luminosities. 

There is some ambiguity about the proper value of (re- 
call L cx f"'"^) to use in blast wave solution for the evolution 
of the quasar light curve, depending on the exact derivation 
from § |2] The solid line in the figure shows our prediction 
for rjL = 18/31, the exact solution to the perturbative accre- 
tion flow in a medium with 7 = 5/3 and kp = (i.e. assuming 
no strong density gradients in the cloud). Allowing for varia- 
tion in 7 and kp in a reasonable range, assuming the black hole 
feedback immediately decouples from the blast wave once the 
blowout begins, or using instead e.g. the estimate of the Bondi 
rate at Rbh in the interior of the blast wave yield iji in the 
range 1/2 < 77^ < 3/2. Therefore, we show in the figure our 
prediction for rji = \/2 (dashed) and = 1 (dotted). 

The primary change is at high accretion rates - a lower 
rjL means a slower decay in the accretion rate, giving more 
time at high accretion rates and high luminosities (see § I4.5> . 
However, the difference in the predictions is small even at the 
highest luminosities, and completely negligible at Mb > -20, 
where the observations are most well-constrained. In fact, the 



systematic uncertainty from different bolometric corrections 
(compare e.g. the open and closed observational points) is a 
larger source of error at all luminosities than this theoretical 
uncertainty in the exact blast wave evolution. 

We also show for comparison the extrap olation of the 
SDSS-2dF quasar luminosity function from iRichards et alJ 
(2005) to z = (see also Bovle et al. 2000), with the Icr range 
indicated as the green shaded range at high luminosities. Sim- 
ilarly, we show the extrapolation of the Ueda et al. ( 2003i) 
hard X-ray luminosity function to z = as open green dia- 
monds over the same luminosity interval, converted to a B- 
band lu minosity function aga in using the bolometric correc- 
tions of iMarconi et al.' ("2004'). We show the corresponding 
(lower right) predictions from th e mode l for merger-induced 
quasar activity of Hopkins et al. (2006a c), which reproduces 
the observed bright quasar luminosity functions over a wide 
range of redshifts and luminosities, as the solid black line 
(adopting the bolometric corrections of IMarconi et alJl2004ft 
and dashed black line (adopting the corrections of lElvisetall 
fT994i) . 

The difference in shape between the quasar and Seyfert lu- 
minosity functions is related to several effects. First, the mass 
function of black holes and, correspondingly, the systems 
"driving" accretion is different (the merger mass function for 
quasars, and the blue galaxy mass function for Seyferts). Sec- 
ond, our feedback-driven model gives slightly different decay 
solutions for both. We explicitly calculate different timescales 
for decay in the two regimes in § |3j and owing to changes 
in the external mass profile and initial energy injection, the 
profile for the decay (i.e. faint-end slope) will be different. 
Essentially, in more violent quasar systems, the feedback 
drives the system to lower luminosities more quickly, result- 
ing in a shallower fain t-end slope of the luminosity function 
( Hop kins et alJ2006cl) . 

The agreement between the predicted and observed lumi- 
nosity functions is good over a wide range of luminosities, 
from Mb > -14 to Mb < -22, and further if we consider 
the predicted quasar contribution from other fueling mecha- 
nisms. At the lowest luminosities Mb > -14, our prediction 
agrees with the [Oil] determinations, but underpredicts the 
luminosity function estimated from Ha measurements. How- 
ever, there are several sources of uncertainty at low luminosi- 
ties. A detailed comparison of the Seyfert 1 and Seyfert 2 Ha, 
[O II], and [O III] luminosity functions in Hao et al. ( 2005), in- 
cluding a comparison between different selection criteria from 
Kewlev et al. (2001) and Kauffmann et al. (2003a), suggests 
that there may be significant contamination by star formation 
at these luminosities in Ha, resulting in a significantly higher 
estimate (note that we compare with the stricter AGN cut the 
authors adopt, although contamination is still possible). 

Additionally, measurement errors and bin-to-bin variation 
are significant at these luminosities, ^ 0.3 dex. Moreover, 
the luminosity function must turn over strongly near these 
luminosities. Although the Seyfert number density implied 
by integrating the observed luminosity function over the ob- 
served range implies that approximately 20% of galaxies host 
a Seyfert in this luminosity interval (an observation necessar- 
ily reproduced by our modeling since our predicted luminos- 
ity and mass functions agree with those observed), extrap- 
olating this luminosity function only 1.5 magnitudes fainter 
would imply that there are more Seyferts than there are galax- 
ies. Finally, alternative fueling mechanisms such as stellar 
winds may become important at these lowest luminosities, a 
point discussed in §|8l 
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Fig. 2. — Predicted intrinsic Seyfert luminosity function from our blast wave model. Solid black line extending to the upper left shows the prediction assuming 
Vl = 18/31 as calculated in §|2| dashed and dotted lines show the predictions varying rii_ t o rii = \/2 and rj^ = 1, respectively. The observed SDSS Seyfert 
luminosity function from narrow-line Ha (i.e. isotropic, intrinsic luminosity function) from lHao et alj RoOSi^ is shown (blue), with both a double-power law 
(squares) and Schec hter function (circles) fi t to the obse rvations, converted to intrinsic B-band using the corrections derived by the authors (filled) and a constant 
conversion followinglMulchae v et alJll994l) :'Elvis et al. ( 1994) (open), each over the observed range. The luminosity function from [Oil] measurements is shown 
in the same style (red). Green diamonds and shaded region show the quasar luminosity function extrapola ted to z = from Ueda et al. 1 2003) and Richards et aQ 
1^05), respectively, with predicted contribution of merger-driven, spheroid-forming quasar activity from lHopkins et aU 1 2006a ) (solid black line in lower right 
for bolometric corrections of Marconi et al. ( 2004), dashed a constant L = 17 Lb). Contribution from stellar winds is shown (dot-dashed) as estimated in §|8] The 
expected Seyfert luminosity function from accretion of cold gas and our model for feedback-driven, self-regulated evolution dominates and agrees well with the 
observed luminosity function from — 15 < Mb < —23, with relatively weak dependence on the theoretical uncertainties for all but the brightest systems. 



4.4. Contribution of Dijferent Morphological Types 

Our analysis enables us to decompose the contributions to 
the AGN luminosity function from different galaxy types. 
At the brightest end, above the break in the extrapolated 
quasar luminosity function, the systems are traditional "bright 
quasars", which may be entering the blowout phase in the fi- 
nal stages of a galaxy merger, expelling the remaining gas in 
a newly-formed elliptical galaxy, a model described in detail 
by Hopkins et al. (2005a-e; 2006a-e). The volume density of 
these objects is very low, <C 10~^Mpc~^, and thus while they 
may dominate the bright AGN population at high redshifts 
they will not be observed locally even in large surveys such as 
the SDSS. 

At luminosities below the break in the extrapolated quasar 
luminosity function and at the brightest end of the AGN lu- 
minosity function, there is a substantial contribution from 
black holes relaxing after the blowout stage in their host 
and black hole-forming mergers. These follow a simi- 
lar decaY_to_die blast \y ave solution described above and 
in iHopkins et al.l J2006bh (both analytically and from sim- 
ulations of mergers) for the specific case of post-merger 
blowout {rii oc f"^), a steeper typical decay than for Seyferts. 
These host galaxies rapidly redden and evidence of distur- 
bance fades quickly, and they will be seen as relatively nor- 



mal ellipticals with large black holes at moderate to low 
accretion rates, with possible evidence for recent (^Gyr) 
merger or star formation activity. This population, with 
precisely these properties and a similar fractional contribu- 
tion to that we predict at the bright end of the local AGN 
luminosity function, is well known observationally (e.g. , 
'Kauffmannet al . "2003 a^. 'Sanchez & Gonzalez-Se rrandl200a 
Sanchez et al. 2004), even in cases where the AGN dominates 
the observed spectrum ( Vanden Berk et al. 2005). 

The bulk of the luminosity range from -15 > Mb > -20 is 
dominated by late-type systems at moderate to high accretion 
rates, fueled by the accretion of cold gas. At the brightest lu- 
minosities, there is some contribution to the AGN luminosity 
function from relaxing ellipticals, but at lower luminosities 
(as is evident from the predicted post-merger AGN luminos- 
ity function prediction in Figure|2j this contribution becomes 
small. Ellipticals do not significantly contribute to the AGN 
luminosity function determined by this fueling mechanism, 
as they do not have a supply of cold gas. The luminosity 
function from this mode of fueling is mainly determined by 
systems of intermediate mass (Mbh ^ 10^ Mq) in Sa/b sys- 
tems and to a slightly lesser extent (owing to their lower gas 
content) by SOs. Sc/d galaxies may not be an insignificant 
contribution to the Seyfert luminosity function, but they do 



Fueling Seyfert Galaxies 



15 



not dominate owing to both their lower characteristic black 
hole masses (by a factor ~ 10-30) and lower number density 
(by a factor ~ 2) compared to both Sa/b and SO systems. In 
gener al, sy stems with small bulges and black holes < 10^ Mq 
(see § 14. 1> may only contribute significant ly at the lowest lu- 
minosities Mb ;J -15. As discussed in § 14.31 above. Sm/lm 
galaxies have such low black hole masses (if any) that they 
contribute negligibly. 

This distribution of AGN hosts for this luminosity inter- 
val is consi stent with observations ( Kauffmann et al.'' 2003at 
ISanchez eT al. 2004; Best et al. 2005), specifically those of 
e.g. Dong & De Robertis (2005) who find that Sa/b and SO 
systems make up most of the low-moderate luminosity con- 
tribution to the Seyfert luminosity function, with a relatively 
small contribution from Sc/d systems. Similarly, the range 
of masses and host galaxy types for which large m > 0.1 
accretion rates are predicted agrees with observational es- 
timates suggesting that present-day black hole growth (in 
the sense of the population of high-Eddington ratio objects) 
is dominated by late-type, seemingly normal systems with 
black hole masses in the range given above (Cowieetal. 
[1996; Steffen et al. 2003; Bai'ser et al. 2003; Ueda et al. 2003; 
iHeckman et al..2004 ). as discussed in more detail in § 14.51 

At the lowest luminosities, there is a substantial con- 
tribution from low-Eddington ratio accretion in late-type 
systems, but also from relaxed ellipticals at low accretion 
rates. Fueling by stellar winds, either from young star 
clusters in systems which still have cold gas or smaller 
contributions from aging stellar populations in old bulges 
has long been recognized as a significant fuel source for ac- 
cretion (e.g., Shull 1983; Mathews 1983; David et al. 1987), 
however the high velocities of these winds yield relatively 
low Bondi rates, and a wide variety of observations further 
show that such systems tend to be accr eting at rates signif- 
icantly below the Bondi estimate ('Fabian & Canizares? 
1^88; Blandford & Beeelman 1999; Di Matteo et al. 
2000; Naravanetal. 2000; Ouataert & Gruzinov 2000; 
DiMatteo et al. 2001: Loewenstein et aL.200L .Bower et aL 
2003r iPellegrini 2005|). 

Nevertheless, these can provide a significant contribution 
to the lowest-luminosity systems, and many "dead" ellipti- 
cals with accretion rates m ^ 10"^- 10"** are expected and ob- 
serve d (e.g.. Ho 2002; Heckman et al..2004;iMarchesini et alJ 
120041 Llester 2005; Pellegrini 2005). These may also be fueled 
by mechanisms other than stellar winds or the accretion of 
hot (virialized) gas, but this seems to account for most of the 
objects here, as calculated in §|8] especially if one accounts 
for feedback r emoving some of th e accreted mass in a steady- 
state solution dSoria et alJ2005bl) . We show the prediction for 
the contribution from stellar wind and hot gas fueling in Fig- 
ure 121 as the dot-dashed line, which demonstrates the lower 
accretion rates of such systems, but the narrower luminos- 
ity range results in a steeper stellar-wind induced luminosity 
function which becomes important at the lowest luminosities. 

4.5. Distribution of Eddington Ratios 

From the evolution of the accretion rate in our blast wave 
solution, we can predict the accretion rate distribution as 
a function of e.g. black hole mass, luminosity, and host 
galaxy properties. Figure |3l shows (upper left) the predicted 
Eddington ratio distribution d//dlogm for "active" (i.e. in 
some stage of blowout) late-type systems determined in Equa- 
tion ( I86> . with the appropriate exponential cutoff at low m 
(such that /d/ = 1). We show this for Mbh = 10^ M©, where 



w 1 (as determined in Equation 1721 ) and typical of the 
black holes which dominate the Seyfert luminosity function, 
but our calculation in § l4.2l demonstrates that this distribution 
depends only weakly on Mbh in late-type systems. 

We show results for three values of rji. rji=\/2 (dotted), 
rjL = 18/31 (solid), and rii = \ (dashed). Because L cx f"'"^, 
larger values of -qi correspond to a more rapid falloff in the 
accretion rate and therefore broader Eddington ratio distribu- 
tions extending to lower accretion rates. In what follows, we 
adopt = 1 8/3 1 , and a different rji will not change the trends 
in Eddington ratio which we find but will systematically shift 
the typical Eddington ratios as shown in the figure. As is clear 
from comparison with Figure|2l the exact choice of rji has lit- 
tle effect on our predicted Seyfert luminosity function, within 
the reasonable range of rji predicted by our blast wave model, 
1/2 < '7l < 3/2. However, rji dominates the systematic un- 
certainty in the estimated Eddington ratio distribution. 

Figure |3l also shows (upper right) the predicted cumulative 
Eddington ratio distribution as a function of host galaxy mor- 
phology. Here, the small differences in the Eddington ratio 
distribution among late-type galaxies are caused by the weak 
dependence of and the duty cycle on host galaxy properties 
(cr and /gas)- The Eddington ratio distribution of ellipticals 
and inactive SOs is estimated from the predicted formation 
times and blast wave decay of merger-induced quasar activity 
(Hopkins et al. 2006a b c). The "cumulative" Eddington ratio 
is, in general, ill-defined, and here we plot the distribution in 
active (i.e. with an event in tn) systems with Mb < -12. 

The lower left panel of Figure |3l shows the (logarithmic) 
mean Eddington ratio and rms dispersion (shaded ranges) as 
a function of Mb from the Eddington ratio PDFs above and 
predicted Seyfert luminosity function in Figure |2] We show 
this for all AGN (black, with yellow shaded range), AGN in 
late-type hosts as modeled herein (blue), and Type 2 AGN in 
late-type hosts as calculated in § 16. II below (red; dispersion 
not shown for clarity but similar to that of all late-type AGN). 
The B-band magnitude Mb represents the intrinsic magnitude 
as in Figure |2] and does not account for extinction. In the 
lower right panel of the Figure, we show the same quantity as 
a function of black hole mass. 

Figure 13 shows our predicted Eddington ratio distribution, 
as in the upper left panel of Figure |3| (but in log-log). Re- 
call that this shows our prediction of d//dlogm, the fraction 
of objects per logarithmic i nterval in rh. NVe compare our re- 
sults with those measured in lYuetalJll2()()l from a sample of 
> 20,000 local SDSS AGN (points), from estimates of [OIII] 
luminosities and black hole masses inferred from measure- 
ments of (T. We show the observations for five values of a, 
from a = 70- llOkms"' (binned by 0.05 dex in log cr). We 
show this range of a because it corresponds approximately to 
Mbh = 2 X 10^-2 x 10^ M©, appropriate for comparison with 
our predictions for Mbh = lO^M©. Furthermore, at larger a, 
there is a significant contribution from "dead" ellipticals (as 
seen in e.g. the black hole mass function of Figure ^ or Ed- 
dington ratio as a function of Mbh in Figure|3jl, and while we 
can predict the combined Eddington ratio distributi on at these 
(T by combining our predictions with those of Hopkin s et alJ 
(2006a b), the relation to our blast wave model and physical 
interpretation is less direct. 

The agreement between our predicted Eddington ratio dis- 
tribution and the observations is good ix^ /v = 0.6 for the 
comparison with the Mbh = 10^ Mq, i.e. a w 100 kms"' data), 
over the observed range. The agreement for 1 /2 < 77/, < 5/6 is 
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Fig. 3. — Upper left: Predicted distribution of E ddington ratios m = M/M^a, for late-type active (i.e. in some stage of "blowout") systems with Mbh = 10^ Mq, 
predicted by our blowout model in Equation 1861 . Results for different rates of accretion rate decay r]L are shown, as labeled (note that the vertical axis is a linear 
scale). Upper right: Cumulative Eddington ratio distribution in active systems (systems having a high-accretion rate event in a Hubble time) with Mb < -12, as a 
function of morphological type, with late type Eddington ratios estimated as above and the E/SO contribution from post-spheroid forming merger decay estimated 
following Hopkins et al. (2006a be). Lower left: Mean Eddington ratio and rms dispersion (shaded range) as a function of Mb from the Eddington ratio PDFs 
and predicted Seyfert luminosity function in Figure|2] for all AGN (black, with yellow dispersion), AGN in late-type hosts (blue; modeled herein), and Type 2 
AGN in late-type hosts (red; dispersion not shown for clarity but similar to that of all late-type AGN). Note that this does not imply higher Eddington ratios for 
Type 2 AGN as a whole, as the luminosity-dependence of the Type 2 fraction is important for such a statistic (see Figurel9l. Dashed line shows the expectation if 
^BH = 3 X lO' Mq= constant (i.e. implying the luminosity function is dominated by differences in Eddington ratio). Lower right: Mean Eddington ratio and rms 
dispersion, in the same style, as a function of Mbh. 



less g ood, but still acce ptable {y^ /v ^ 1 .5-2.0). Note that we 
follow !^ et al.l 12005") and exclude the lowest Eddington ratio 
bin for each a owing to significant incompleteness. In order 
to compare with their observations, we have convolved our 
prediction with the distribution of [OIII] luminosities (i.e. al- 
lowed for a small dispersion in bolometric correction), which 
is why our prediction is non-zero (although small) above an 
Eddington ratio of unity. The predicted Eddington ratio distri- 
bution in this mass range depends only weakly on ct, evident 
in e.g. our Equation \12\ and in Figure |3l and therefore the 
observations over this entire range in <j all agree with our pre- 
diction in a nearly self-similar manner. 

That the observed Yu et al. (2005) in distributions agree 
with our exact blast wave solution with rji = 18/31 and rea- 
sonably well with e.g. the range of solutions 1/2 < rji < 5 /6 
is somewhat surprising. The slope this implies (i.e. loga- 
rithmic slope of d«/dlogm) of -l/ry/, = -31/18 is substan- 
tially different from the slope of « -0.8 fitted to the low- 
m behavior of the m distributions by the authors, which 
would imply ?/£ sa 1.26 ± 0.1. The authors use this fit to ar- 
gue for a self-similar quiescent disk evolution model to ex- 
plain the late-time L cx f"'"^ behavior However, this fit from 
|Yu et al. (2005) is derived primarily from the low-;7i behavior 
of their a « 200km s"' objects, which we have noted (fol- 



lowing numerous observational studies) are a different pop- 
ulation, dominated by ellipticals decaying from larger accre- 
tion rates or fueled by e.g. stellar winds and diffuse hot gas 
(e.g., Kauffmann et al. 2003a; Pellegrini 2005). Furthermore, 
the authors acknowledge that this does not provide a good 
fit to objects with larger a, which is surprising since these 
objects are otherwise observed to be a similar AGN popula- 
tion to those with a w 200km s"' (i.e. both are dominated by 
early-type systems). The authors also acknowledge that this 
model for the activity predicts that the active black holes will 
be observed at masses substantially below their final (Mbh-c 
relation) masses, by about Idex, whereas the observations 
(e.g., Barth et al. 2005) suggest the opposite trend. In §|5lwe 
demonstrate that our model predicts this observed tendency 
rather than smaller black hole mass es. 

Additionally, the model fitted bv lYu et all J2005h does not 
provide a good match to the high-;« distributions, even for 
those values of cr for which the overall fit is most acceptable, 
giving e.g. typical x^/'^ ^ 4-6 for a ^ 100-200 kms"' (at 
m > 10"^). Over the entire observed range, our model pro- 
vides an improv ed fit, x^/i^ » .6 compared to x^/^ w 1.6 
from the fits of lYuetalJ (^2005); alternatively, Ax^ w 11.0 
for e.g. (T = llOkms"^ These arguments suggest t hat while 
the self-similar quiescent disk evolution discussed in lYu et all 
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Fig. 4. — Predicted Eddington ratio distribution of active Seyferts from our blast wave model in §|2|and for estimated duty cycles in $ 14.21 as in Figure|3](upper 
left), for our exact solution rjt = 18/31 and the range jjl = 1 /2, 1 as labeled, for black holes with Mbh = IO'Mq (although this depends only weakly on Mbh)- 
The observed accretion rate distribution from Yu et al. ( 2005) is plotted for each of several values of a as labeled, con'esponding to Mbh ~2x 10*— 2x lO' Mq . 
Both the power law-trend at high m (dependent on the feedback-regulated light curve decay) and flattening/turnover at low m (dependent on excitation rates and 
duty cycles) in the observed distribution agree very well with the predictions of our model (y^ jv = 0.6), whereas models without feedback-driven self regulation 
of the accretion rate produce much flatter distributions. 



( I2005 [) may be important for some range of early-type sys- 
tems at low accretion rates, the observations in the range of 
black hole masses shown in Figure |4] are better explained by 
our model for feedback-driven blast wave evolution. 

The turnover at low-m in the Eddington ratio distribution 
is, for any power-law distribution of accretion rates, set en- 
tirely by the duty cycle at high m. In our language, for a given 
slope rji, this turnover occurs at a value in^i„ determined by 
Equation (I88> as a function of the duty cycle (i.e. lifetime and 
rate of activation) at high m. There is some uncertainty in the 
exact normalization of this turnover, as our theoretical predic- 
tion depends on scaling arguments and the detailed numerical 
prefactors are sensitive to the local kinematics near the black 
hole. Therefore, our prediction for the location of the turnover 
is more uncertain than our prediction of rji. 

Currently, the observations do not resolve the turnover. Do- 
ing so would constrain the duty cycle as a function of m; for 
example, the difference between a sharp turnover as we have 
assumed and a shallow (e.g. symmetric power-law) turnover 
for the »Jn,in of our r//, = 18/31 prediction amounts to a fac- 
tor of ~ 2 change in the duty cycle Sq. Therefore, observa- 
tions of the shape of the accretion rate distribution a t a given 
a, probing only moderately fainter luminosities than lYu et al.l 
(|2P05), can constrain this quantity without the large system- 
atic uncertainties inherent in determining the normalization of 
the distribution when the turnover is not resolved. We note, 
however, that at least the flattening of this distribution at low 



m is resolved by the observations of lYu et all J2005I) - if our 
prediction were simplified to a pure power-law, it would over- 
predict by ^3-6(7 each of the lowest four m observations 
plotted. 

Although incompleteness in the lowest-m bin at each a may 
be significant, the observations appear to favor the turnover in 
the Eddington ratio distribution moving to larger m at lower cr, 
an expectation broadly consistent with our result in the lower- 
left panel of Figure |3| Figure |5] shows our prediction for the 
mean Eddington ratio as a function of black hole mass from 
Figure|3] compared to the estimates of Heckman et al. (20041) 
(squares). Neither the observations of Heck man et alJ (120041) 
nor those o f f^ et alJ t2005) resolve the turnover in accretion 
rate distributions at low-m, and therefore cannot determine a 
"mean" accretion rate in a proper sense. Therefore, we in- 
stead consider the cumulative distribution of accretion rates 
from the observations for Mbh ~ 3 x lO^Msun, and fit this to 
a distribution of Eddington ratios similar to our prediction in 
Figure 13 above, but allowing both the slope and duty cycle to 
vary freely. From these fits, then, we can constrain the peak 
and mean of the Eddington ratio distribution. The vertical er- 
ror bars of the point at Mbh = 3 x lO^Msun in Figure |5] show 
the 1(7 allowed range from these fits, and reflect the system- 
atic uncertainty in the absolute value of the mean accretion 
rates at a given Mbh or cr determined observationally. 

While the mean Eddington ratios are not well-determined 
observationally, the trend of Eddington ratio with mass is 
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log,o(MBH/MQ) Eddington Ratio 

Fig. 5. — Left: Mean Eddington ratio as a function of black liole mass from Figure |3] compared to estimates from tlie observations of 'Heck man et alj 
120041) (squares). The point at Mbh ~ 3 X lO^Msun sh ows the absolute systematic uncertainty in the mean Eddington ratio distribution from fitting to the (high 
Eddington ratio only) distributions of 'Heckman et alj 12004) . subsequent points show the relative uncertainty in the trend of Eddington ratio vs. Mbh. The 
observed trend of decreasing Eddington ratios with Mbh is reproduced, as larger-MBH populations are increasingly dominated by large spheroids without cold 
gas supplies for stochastic accretion. Right: Tentative Eddington ratio distributio n from Figure|3](upper left), summed over all morphological types, compared to 
the observationally estimated Eddington ratio distribution from Marchesini et al. 120041) (histogram). Note, however, that the vertical axis here is linearly scaled, 
and the predicted bimodality can easily be erased by e.g. different values of til. 



more well-constrained. Using the best-fit absolute normal- 
ization from the fit to the Mbh = 3 x lO^Msu,, data, we show 
the expected mean Eddington ratio as a function of Mbh at 
Mbh = 10^ 3 X 10^ 10^ 3 x 10^ and 10'^ M© from the ob- 
served trends in Heckman et al. ("2004") (see their Figure 3). 
Here, the vertical error bars represent the relative error in Ed- 
dington ratio as a function of mass. Thus, the typical mean 
Eddington ratios do agree with our predictions, but the un- 
certainties are large because the downturn of the m distri- 
bution is not measured observationally. However, the trend 
of mean Eddington ratio with black hole mass is reasonably 
well-constrained, and agrees with our predictions for black 
hole masses Mbh = 3 x 10*'- 10'' Mq. 

The trend arises because for higher black hole masses, the 
population is increasingly dominated by ellipticals which do 
not have a supply of cold gas. These systems are presumably 
either fading from earlier bright quasar activity or accreting 
in a quasi-steady state from virialized hot gas or stellar mass 
loss, at accretion rates well below those of late-type Seyferts. 
Furthermore, the mean age of ellipticals may increase with 
mass, here calculated following the determination as a func- 
tion of bulge and black hole mass in Hopkins et al. (2006b), 
further contributing to the trend of decreasing m with increas- 
ing Mbh even in pure elliptical systems. 

From Figure |5] our modeling predicts that most of the 
black hole mass growth in the local Universe occurs in 
relatively low -Mbh syst ems with Mrh 5. 10 ^ M^, consis- 
tent with observations ( Heckman et al.l I fool . Of course, 
observations of quasars, compared to our predictions for 
the e volu tion of cloud-fueled activity in late-type galaxies 
in § 14.61 suggest that this is not true at higher redshifts, 
where more massive black holes are formed. Although 
this c onstitutes the low-r edshift end of "cosmic downsizing" 
(e.g., iCowie etaL 119961) seen observationally in AGN and 
quasar evolution (iPage etaPl 19971 iMivaii et alJl2000l 1200 U 
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l2005HHasinger. Mivaii. & Schmidi2005l) . our modeling does 
not imply that the local activity is fundamentally cosmolog- 
ical. The transition to lower black hole masses growing at 
lower redshifts through black hole growth via merger-driven 
quasar fueling is cosm ological and traces galaxy downsizing 
("Hopk ins et alJl2006b ). But, the relatively high rates of ac- 
cretion in late-type galaxies compared to higher-MBH systems 
predicted here is more a function of the high-MBH ellipticals 
fading with time, rather than any sudden "turning on" of the 
Seyfert population. 

In Figure |5] we also compare our result for the total distri- 
bution of accretion rates across all morphological types from 
the upper right panel of Figure |3] (solid line) to that deter- 
mined by Marchesin i et al, ( 2 0(14,) (shaded hi st ogram) , which 
agrees with the estini ates of e.g. |h3 ll2002t : iTesteri ^2005): 
'Hopkins et al.' ('2005e') and extends to lower m than that of 
e.g. Yu et al. (2005) and includes inactive systems. The pre- 
dicted Eddington ratio distribution is bimodal, with Seyferts 
comprising the majority of the high-m peak and relaxed ellip- 
ticals dominating at low m ^ 1Q~* accretion rates. There is 
broad agreement, and the slight horizontal offset between the 
theory and measurements is within the systematic uncertainty 
of the observational estimate of m. 

While this agreement is suggestive, we caution that the 
predicted bimodal distribution is not firm. The vertical axis 
shows a linear scale, and, from Figure^] a small change in our 
choice of rji or the duty cycle at high m (which determines the 
'«min where the high-;?! distribution turns over) would broaden 
the distribution to < 10"^, erasing the bimodality. Further- 
more, our prediction depends on the Eddington ratio distri- 
bution of quiescent ellipticals, estimated from their accre- 
tion rate decay according to our blast wave model following 
mergers and simulations of merger-driven black hole growth 
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iHopkins et alJl2006ald) . but once these low accretion rates 
are attained, other processes may dominate the gas inflows 
and accretion, and thus our estimate may not be appropriate 
for this end of the distribution. 

In the lower-left panel of Figure |3j we show our predic- 
tion (solid line) for the mean accretion rate as a function of 
luminosity, compared (dashed line) to the corresponding Ed- 
dington ratio for a black hole of mass 3 x lO'' Mq at each Mb 
(i.e. approximately a logarithmic slope of -1/2.5). This slope 
is what would be expected if the observed luminosity function 
were purely an Eddington ratio sequence. Although the dif- 
ferences between our predicted distribution and that shown 
are important, the rough similarity between the two, at least 
for Mb > -20, implies that the range of the observed luminos- 
ity function is dominated by differences in Eddington ratio, 
and not by different black hole masses at similar relative ac- 
cretion rates. Since e.g. host galaxy mass and luminosity are 
correlated with black hole mass, we expect that these quan- 
tities will not be strongly correlated with the observed AGN 
luminosity, as a wide range of AGN luminosities can have a 
similar distribution in constituent black hole masses. This has 
been seen in observ ations of local AGN and quasars and their 
host galaxies (e.g.. 'Bahc all et aljn997: .McL ure et al .lJ1999|: 
[Hamilton et al. 2002; Wo o & Urrvi2002tlO'Do wd et a l.l20oi 
pahnke & Wisotzkiii2003tlHao et alJl2005HVanden Berk et alJ 

Note that the prediction for the Type 2 Eddington ratios in 
this panel, while slightly above those of Type 1 objects at 
high luminosities, does not imply that we predict typically 
high er ac cretion rates in Type 2 Seyferts. Rather, as shown 
in § 16.11 our modeling indicates that Type 2 Seyferts are a 
small fraction of the population at these high luminosities, but 
a large (comparable) fraction at low luminosities. Therefore, 
our modeling predicts that the characteristic accretion rates 
of Seyfert 2s should be lower than those of Seyfert Is. In 
general, the differences in Eddington ratio distributions be- 
tween Type 1 and Type 2 systems is a more sensitive function 
of their characteristic luminosities and determined by e.g. the 
dependence of the Type 2 fraction on luminosity, than it is a 
difference between the two populations at a given luminosity. 
(This refers to intrinsic luminosity, as Type 2 systems will, by 
definition, be extinguished at some frequencies.) 

4.6. Redshift Evolution of Seyfert Activity 

Our estimate of the duty cycle as a function of m depends 
primarily on the gas fraction and cr, and the convolution to 
give a luminosity function depends on the black hole mass 
function in late-type galaxies. From the evolution of these 
quantities it is straightforward to determine the evolution in 
this luminosity function. 

Because the mass gain ed in a "blowout" event is small (at 
most a factor ^ 2; see § \5A\ . we determine the black hole 
mass function at each redshift of interest from the evolution 
of the late-type galaxy mass function, using the method of 
§ 14.31 From the spectral-type separated mass functions of 
IFontana et al. ( 2004) determined up to z = 2 from the K20 sur- 
vey, we expect little evolution in the late-type mass function 
up to z ^ 1, and then a decrease in the number density of these 
galaxies. In detail, we trace this evolution up t o z = 0.5 using 
the K and /-band luminosity functions of e.g. 'F eulner et alJ 
(|2003); Pozzetti et al. (2003); Dahlen et al. (2005). Normahz- 
ing the mass functions to be the same in the overlapping red- 
shift interval 0.2 < z < 0.55, we extend this with the morpho- 
logically classified spiral mass functions from iBundv et alJ 



( I2005h and lFranc eschini et al l (l2006h up to z = 1.4, and then 
repeat this with the Fontana et alj (120041) observations to infer 
the mass functions to z = 2. 

These estimates are roughly consistent with the B-band evo- 
lution in late-type luminosity functi ons which is more well- 
const rained (e.g., Id e Lapparent et al. 2004; Giallong o et alJ 
l2005t iFaber et al. 2005), measured evolution in the c umula- 
tive m ass function ( de Propris et aL.1999i: .Drorv et alJl2004 
and evolution in black hole mass functions from 
spheroid luminosity functions and integration of the continu- 
ity equations for AG N luminosity functions ( Yu & Lu 2004i 
Tamura et all 120051) . Given this, our approach appears to 
be reasonable, but the observations have large uncertainties. 
However, we find that there is little evolution in the Seyfert 
luminosity function in any case, and these uncertainties are 
not important. 

We estimate the evolution in the typical gas fractions of 
disks by assuming that gas is consumed on a timescale 
related to star formation fsF ^4-8 Gyr, giving /gas ^ 
exp(-f//(z)/fsF), where tniz) is the Hubble time at red- 
shift z. We take fsp ~ 6 Gyr, as this gives a Milky- Way 
like /gas ~ 0.1 at z = 0. This prescription a nd charac- 
teristi c timescale follow from observations (e.g., 'Kennicutj 
119981 Rownd& Young 1999; Martin & Kennicutt 2001) 
cosmological simulations (e.g., Springel & Hernauist 200^ 
[Hernquist & Springel 2003; Nagamine et al. 2004, 2005a b), 
simulatio ns of star forrnat ion rate evolution in isolated disk 
galaxies dLi et alJl20053g) . and comparison of these simula- 
tions and the predictions for quas ar activity with merger rates 
and merger luminosity functions ("Hopk ins et al.l l2006e). As 
gas fractions increase, the event rate for cold gas accretion 
can increase, producing a higher duty cycle at large m from 
Equation i&ll . Although our choice for the evolution of the 
gas fraction is somewhat arbitrary, we find a sufficiently weak 
dependence that observational uncertainty is not significant. 

Figure |6l shows our predicted Seyfert luminosity function 
(black lines with no shaded range) at various redshifts. The 
solid lines show the full prediction, accounting for both evo- 
lution in the late-type galaxy mass functions and the evolu- 
tion in typical gas fractions of disk galaxies. The dashed line 
shows an estimate of "maximal evolution," i.e. allowing the 
gas fraction to increase with redshift, but not allowing any de- 
crease in the mass functions. The dotted line shows the z = 
prediction from Figure|2lat each redshift - i.e., no evolution. 

For comparison, at each redshift, we show the fitted quasar 
luminosity function of Ueda et al. (2003) (black solid line 
with yellow range showing the uncertainty owing to different 
bolometric corrections), converted to Mb as in Figure El The 
"low" (0.2 < z < 1) redshift behavior of this luminosity func- 
tion has also been compared with the luminosity functions 
of e.g . Richards et al. (2005); Hasinger. Mivaii. & SchmidJ 
(^2005) which have more objects at these redshifts, and they 
give consistent results. Alternatively, we could plot the expec- 
tation for merger-driven quasar activity calculated from the 
models and simulations of quasars in mergers from Hopkins 
et al. (2005a-e; 2006a-e), which is essentially identical to the 
plotted quasar luminosity functions. 

From Figure |6l we see that the quasar luminosity function 
evolves rapidly, while our predicted Seyfert luminosity func- 
tion hardly evolves even out to z = 2. While our estimates 
of the evolution of e.g. gas fractions and late-type mass func- 
tions are uncertain, the similarity (relative to the evolution in 
the quasar luminosity function) of the expected evolution, no 
evolution, and maximal evolution cases means that unless an 
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Fig. 6. — Predicted Seyfert luminosity function owing to the stocliastic fueling, at several redshifts (as labeled). The observed quasar luminosity func- 
tion from Ueda et al. (2003) (black line with yellow shaded range) is shown, converted to Mb as in Figure l2l The solid line ( no sha ded range) shows 
our prediction for the Seyfert luminosity function, determined from the evolution in late-type mass functions ide Propris et al. 1999 HDrorv et al . 2004, 2003 
^ntana et al. 2004; de Lapparent et al. 2004; Bundv et al. 2005; Dahlen et al. 2005; Giallongo et al. 2005; Franceschini et al. 2006) and gas fractions i KennicutJ 
^98; Rownd & Yo unail999i;iMartin & Kennicutliii200kiSDri ngel & Hernauist 2003; Li et alj.2005^h.) . Dashed lines shows the "maximal evolution" (i.e. al- 
lowing gas fractions to increase with z but not allowing any mass function decrease). Dotted lines show the z = prediction of Figure |2] i.e. no evolution. 
Regardless of our assumptions, the evolution in the luminosity function of systems fueled by quiescent, non-merger driven cold gas accretion is weak, and the 
quasar luminosity function at the redshifts of peak quasar activity (and quasar luminosity density) must be dominated by other fueling mechanisms. 



extreme model for late-type galaxy evolution were adopted (at 
odds with the observations), there would be essentially no dif- 
ference in the predicted evolution of this luminosity function. 
Thus, these observational errors are not a significant source of 
theoretical uncertainty for the key conclusions from Figure|5] 

From the figure, it appears that our prediction may be sys- 
tematically high, by about a factor of ~ 2, at the lowest lumi- 
nosities, which is possible, given the uncertainties in our esti- 
mate of the duty cycle. If the normalization at these luminosi- 
ties is matched to the extrapolated faint end of the quasar lu- 
minosity function, then it implies an even smaller contribution 
from this activity at these redshifts. However, it is important 
to note that the faint end quasar luminosity function plotted 
here is extrapolated below the observations at the higher red- 
shifts, and may be incomplete at the lowest Mb ^ -18 lumi- 
nosities, corresponding to or slightly belo w the faintes t bins i n 
even the deep luminosity functions of e.g. 'Ueda et al.'('2003"); 
iHasinger. Miyaji. & Schmidt (2005); La Franca et al. (2005) 
at low (z < 1) redshifts. 

For local, z = AGN, the contribution from stochastic 
cold gas accretion dominates the luminosity function up to 
Mb ~ -22, near the traditional division between Seyfert and 
quasar activity (see Figure |2j- However, above z 0.1, this 
shifts to lower luminosity Mb ~ -20, and the relative contribu- 
tion at higher luminosities steadily decreases. This luminos- 
ity function does still represent a significant contribution to 



the faint end (i.e. below the "break") of the quasar luminosity 
function up to z ~ 0.5. However, by z = 1, the "Seyfert" (i.e. 
our predicted) contribution to the quasar luminosity function 
is only significant at about 4-5 magnitudes (i.e. two orders 
of magnitude in luminosity) below the break, and by z = 2 
only important about 7 magnitudes, or 3 orders of magnitude 
in luminosity, below the break. Thus, at the epochs gener- 
ally associated with "quasar activity," the contribution from 
this accretion mechanism is negligible, even at luminosities 
far below the deepest observations. 

That the quasar luminosity function evolves more rapidly 
than our prediction for the Seyfert luminosity function indi- 
cates that different modes of fueling are likely responsible for 
the two populations. Indeed, merger-induced quasar activity 
can account for the evolution in the quasar luminosity func- 
tion, and merger rates of gas-rich galaxies evolve more rapidly 
than the mass function of gas-rich galaxies. Note that at any 
time the number of gas-rich galaxies undergoing mergers is 
small, so a relatively large number of mergers providing the 
observed quasar activity does not demand strong evolution 
in the disk galaxy mass function. Hopkins et al. (2006b ^ 
study the evolution of quasar, red galaxy, and merger lumi- 
nosity functions and argue that a self-consistent mapping be- 
tween the two implies that a significant fraction of the bright- 
est quasar activity must result from mergers. 

We can use our predictions to determine the contribution of 
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quiescent fueling to the buildup of the black hole mass den- 
sity of the Universe and cosmological backgrounds, such as 
the cosmic X-ray background. From the analysis in § |5l or 
the mean Eddington ratio plotted as a function of Mbh in Fig- 
ure |3l this process will contribute negligibly to the growth 
of black hole masses above a ~ lOOkms"'. Larger-mass 
black holes ( Mrh lO^Mfr^) dominate the present black hole 
mass density ( Salucci et al.''1999'; 'Marconi & Salvati"2002^ 
IYu & Tremaine 2002; Ferrarese 2002; Aller & Richstone 
|2p02; Shankaret al. 2004), and thus the contribution of qui- 
escent fueling to the total black hole mass density is negligi- 
ble. Even at relatively low Mbh = lO*"- IO^Mq, the quiescent 
mode of fueling will increase black holes masses by only up 
to a factor ~ 1.2-2.0 (see § 15. 4> . This can also be seen in 
Figure n where the dashed line shows the inferred BHMF 
if we allow for the modified Mbh - c relation predicted from 
this fueling mechanism in § |5j in other words allowing for 
this fueling mechanism to grow small-mass black holes sub- 
stantially. Because this also results in an increased scatter in 
the relation, the effect is "blurred out," and we see in Fig- 
ure n that the net effect on the black hole mass function is 
small, within observational uncertainty. Therefore, while the 
contribution to the mass growth of low-mass black holes is 
non-negligible (although probably insufficient to grow them 
by orders of magnitude from much smaller seeds), the contri- 
bution to the integrated black hole mass buildup and cumula- 
tive black ho le mass density is small, satisfying e.g. the con- 
straints from i Hopkins et al] ( l2006 a) that non-merger driven 
fueling mechanisms are relatively unimportant in these inte- 
grated quantities. 

Given this small contribution to the integrated black hole 
mass density, it follows from energetics that the contribu- 
tion from this fueling mechanism to the X-ray bac kground is 
also small. Although observations (e .g., iMivaii e t al. 2003 
lUeda et alJ l2003t iBarger et alJ 12005ft have suggested that 
much of the X-ray background is built up at z < 1, where this 
fueling mechanism is relatively more important, we find a sig- 
nificant contribution only at faint luminosities, which only ap- 
proach the break in the quasar luminosity function for z ^ 0.5. 
This is consistent with our previous estimates for the X-ray 
background spectrum determined by merger-induced quasar 
activity in Hop kins et al.i (2006a), which also suggests that 
other mechanisms should not dominate the X-ray background, 
even at redshifts z < 1. This distinction cautions against ex- 
trapolating models for obscuration based on local AGN (e.g. 
traditio nal toroidal models) in X-ray background synthesis 
(e.g., Setti & Woltier'1989':'Madau et al.'1994VComast ri et"al] 
fl995.:.Gini et al. 1999, 2001; Treister & Urry 2005), as obscu- 
ration mechanisms in these different populations are distinct 
physically ( Hop kins et aLi2005bid) . and indeed are observed 
to follow quantitatively different, albeit qualitatively s imila r, 
trends in obscuration as a function of luminosity (see § I6.1> . 

5. BLACK HOLE GROWTH AND THE Mbh - ^ RELATION 

5. L Slope of the Mbh ~ f Relation 

In our scenario, when a molecular cloud of mass en- 
counters a black hole of mass Mbh, there will be a short 
period of high accretion. For a cloud moving at a speed 
~ cf^"^ ~ lOkms"', the Bondi accretion rate will be 

. A-KaG^Mmp^\ts ^ ^ ,„^, 

and accretion will be Eddington-limited with in= 1 . The ac- 
cretion rate will drop rapidly once feedback energy unbinds 



the cloud and heats or expels the surrounding gas. 

If feedback impacts the gas energetically, then the 
"blowout" phase begins roughly when the coupled radiant en- 
ergy is comparable to the cloud binding energy. The system 
will follow the "blast wave" solution if it accretes at rii « 1 for 
a time Af ^ f„-, ^ Rc\/<J (see § I3.2> . and an energy sufficient 
to unbind the cloud (i.e. a fraction ^ 1 of the binding en- 
ergy of the cloud) couples to the gas in a shorter time. This 
timescale is less than or comparable to the dynamical time of 
the cloud ^ 1 1 \fGp and the cloud crossing time Rjd^^^^ (i.e. 
timescale of interaction between the black hole and bulk of 
the cloud), so it does not matter which we use in our analysis. 
They are also <C fs, the Salpeter time, so we are justified in 
assuming that the black hole mass is approximately constant 
and the accretion rate is Mfidd = Mbh/Js over this interval. The 
blowout criterion then becomes 



Mbh 9 ~ 9 

r]LAt = 'qe,- tri,c^>fh Md(t)<J 



(92) 



where rj is the feedback coupling efficiency and is a numer- 
ical coefficient which depends on the bulge profile and gives 
the binding energy at r = as a function of the observed a (i.e. 
line-of-sight averaged velocity dispersion within the effective 
radius). For a Hernauist ( 199QJ spheroid, (f) = 10.1, so this 
factor is significant. This then determines a minimum black 
hole mass for blowout of 
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(93) 

-lOOcm-V VlOOpcV \7]T^/fi,J 

where r* is defined in Equation J72> . If the black hole is ini- 
tially formed along with the bulge as is expected from models 
of the Mbh - cr relation with Mbh c>c cr"*, it may therefore have 
to grow by a non-negligible factor at small-cr before a blowout 
occurs, and this implies a slightly shallower Mbh - cr relation 

at low masses with Mbh oc a . 

I iu' above derivation is similar to that of ISilk & Reej 
11998), except that the characteristic timescale of an event 
and total gas mass are set externally and depend differently 
on (7, from e.g. the dynamical time and gas mass in a merger. 
However, it may not be the case that an energy criterion deter- 
mines the Mbh - o" relation. For "galaxy-scale" co-formation 
of bulges and black holes in mergers, the bulge mass goes as 
M cx cr'*, and the dynamical time goes as oc a/cr cx cr, suggest- 
ing Mbh oc (t^, so the actual relationship may be shallower. 
A momentum-based coupling and "blowout" sug- 
in better agreement with observations 
For momentum coupling, the injected 
momentum is ryi/c, and the criterion becomes 



gests Mrh oc & 



i A Mbh 

77- Af = TItr fm C > Md tJ , 
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which gives a minimum black hole mass 



(94) 



Mbh = 1.13a X lO'^Mg 



(T 



lOOkms-1 
2 / 1.0 



" (lOOcm-0 (lOOpc) (77T*//J' ^^^^ 

where the more efficient coupling (77 ~ 1) for momentum - 
driven feedback is adopted following Murray et al. ( 2005!), 
when the molecular gas is optically thick to the AGN radi- 
ation, at least outside the dust sublimation radius. 
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If the typical spheroid M- cr relation is given by 

M\,y. = 1 .35 X IO^Mq f -\ (96) 

then the minimum black hole mass required to expel a molec- 
ular cloud is greater than Mgj^ at a mass of ~ \^ Mq, or 
G w lOOkms"'. In general, the Mbh-ct relation at any cr 
will be given by the larger of either the "bulge-formation" 
Mbh oc ct"^ or the relation determined above, Mbh oc (or 
Mbh oc cP' for momentum coupling). Thus, at sufficiently low 
masses, the shallower relation predicted above will dominate, 
yielding a shallower slope to the Mbh - c relationship and in- 
troducing a break or some curvature in the relation. 

The location of this break is only loosely constrained since 
the coupling efficiency is not well-determined, and we use 
the normalizations above as they give a good match to ob- 
servations for the merger-induced Mbh oc cr'* relation (e.g., 
|Di Matteo et al. 2005; Murrav et al. 2005). However, the 
slope and scatter (see § 15.21 below) induced by this fueling 
mechanism are robust predictions independent of these uncer- 
tainties, and the "break" location is reasonably constrained to 
^ 10^- IO^Mq based on estimates of the feedback efficiency 
from t he normalization of t he merger-driven Mbh oc a'^ rela- 
tion ( Pi Matteo et alJl2005 l). although this does not apply if 
e.g. the accretion mode and coupling change at lower masses 
or luminosities. 

In Figure we show (solid line) the predicted relation 
for both merger-drive n co-formation of bulges and black 
holes with Mbh oc a"^ (iDi Matteo et al.''2005). matchi ng the 
observed relation (jp errarese & Merritt 2000; Granat o et alJ 
12000- "Tre maine et al . 2002), and for our Seyfert molecular- 
cloud mode of accretion as predicted by Equation ( I93> for 
energy-driven coupling (dashed line) and Equation ( I95> for 
a momentum coupling (dotted line). For comparison, ob- 
servations are shown from various compilations as indicated 
in the figure caption. Note that the (highly uncertain) re- 
verberation region geometric normalization factor adopted in 
iQreene & Ho ( 2005a) dif fers shghtly (by ^ .2 dex) from that 
of >Onken et al. (.2004); .Nelson et all j2004 ; iPeterson etTD 
j2004l) : lRarthetal.l (120051) ." Since we cannot predict the nor- 
malization of the Mbh - c relation to this accuracy and are, 
in any case, interested only in the slope and scatter of the re- 
lation, we adopt the same normalization for all the samples 
used. 

The comparison between the predicted energy-based rela- 
tion and the low-mass observations is suggestive. We con- 
sider the 37 ob served sy stems below a » 100kms ~\ pri- 
mai-ily from the Earth et al. (2005) and Greene & Ho (2005a) 
samples but also including Pox 52, NGC4395 (for which we 
adopt the upper limit a = 30kms"', as a lower value will only 
further favor our estimated Mbh - cr relation over the extrap- 
olated Mbh oc a * relation), and fou r [jremaine et al. ( 2002) 
objects and four Onken et al. (2004) objects. For these ob- 
jects, if we ign ore intrinsic scatter in the relationship, the re- 
duced of the lTremaine et alJll2002^ fitted relation is x^/j^ = 
16.8, whereas our prediction gives x^/'^ = 12.1 (an abso- 
lute Ax^ ~ 174). These are both unacceptable fits, however, 
so allowing for a constant intrinsic scatter o f 0.27 dex gives 
X^lv= 1.72 (marginally acceptable) for the iTremaine et all 
(|2002) fit, compared to x'/i/ = 1 .02 for our prediction (abso- 
lute Ax^ w 26, still highly significant). Further allo wing for 
the increased intrinsic scatter at low a we predict in § l5.2l be- 
low, this drops to x^ jv = 0.21. Above this cr, the comparison 



is not meaningful because we do not predict any substantial 
correction to the Mbh - cr relation. 

The observations at low black hole masses appear to fa- 
vor our prediction over a pure bulge and black hole co- 
formation (i.e. merger-driven) Mbh - cr relation, which yields 
a scahng Mrh oc cr'*, roughly i ndepende nt of mass (see e.g., 
iDi Matteo et ialJ2 005; Roberts on et al.l200 5). However, there 
is only a marginal (~ 1 - 2cr) detection of this trend, and both 
the measurement errors and systema t ic uncertainties (e.g. , 
IPeterson eta r 2005; Earth et al. 2005; Greene & Ho • 2005bh 
are large at low masses. Future observations of the Mbh-ct 
relation at low masses will provide a strong test of our model, 
even potentially distinguishing momentum coupling vs. en- 
ergy coupling models. 

5.2. Scatter in the Mbh ~ cr Relation 

If molecular cloud accretion is important for the growth 
of low-mass black holes, this mode can influence the scat- 
ter in the observed Mbh - cr relation. For example, the crit- 
ical Mbh for entering the "blowout" phase in Equation \95\ 
is proportional to /?^[ cx M^l^ , so a spectrum of sizes of 
molecular clouds will introduce scatter in the Mbh - cr re- 
lation. This scatter will be significant if it is larger than 
the ^ 0.3 dex scatter intrinsic to the Mbh-ct relation from 
spheroid-black hole co-formation. The observed scatter in the 
relationship is approximately constant at high masses (e.g., 
5'remaine et al. 2002; Novak et al. 2005), and simulations of 
merger-induced spheroid formation also predict the scatter to 
remain roughly constant at low masses (.Di Matteo et al. 200^ 
l^bertson et al. 2005|^ 

As discussed by ^McKee & Ostrikei' ('1977'), observations 
of co l umn density distributions in the galaxy (e.g., Hobb3 
1974"; 'Payne et al. 'l983; Welty et^ 11994 ). LMC (e.g., 
Oestreicher & Schmidt-Kaler 1996) and AGN hosts (e.g., 
Tytler 1987; Hopkins et al. 2004) indicate that the number of 
clouds per unit column density varies as approximately N^, 
implying (for a constant cloud density) that the number of 
clouds with radii in the range R ^ R + dR per unit volume is 
cx R~* dR. The rate at which clouds collide with the black hole 
is 

I^Wevent ^ fi^^jj-Rl, Vd OC — dluR OC — \-prdlnMcl, (97) 

dt R mJ 

where n^i oc R~'^ dR and the independence of Vd and R are used 
in the second equality. Therefore, the +lcr (68% R < R') colli- 
sion probability corresponds to a cloud about four times larger 
in radius than our "typical" cloud, i.e. M^i ~ 70 times bigger 

2/3 

than normal, orM^j ^ 17.1 times normal. The +lcr disper- 
sion in the Mbh - cr relation should then be w 1.2 dex at the 
lowest Mbh, larger than the ~ 0.3 dex intrinsic to the Mbh - cr 
relation determined by co-formation of bulge and black hole. 

The lower-MBH limit to the molecular cloud accretion 
Mbh - cr relation will be determined by extrapolating the 
Mbh - cr relation for co-formation of spheroids and bulges. 
Although technically the PDF for encountering a cloud of 
size R (Equation 1971 ) can be estimated down to the minimum 
cloud size ^ 0.5 pc, this is unimportant for the scatter in the 
A^BH - cr relation. For/?<C 10-lOOpc (i.e. less than "typical" 
large cloud sizes), the minimum black hole mass to "blowout" 
from Equation (|95} is smafl, e.g. Mbh ^ 10"* M© for/? ^ lOpc 
(and recall that we are considering intervals in \ogR). Thus 
the contribution to the Mbh - cr relation from these small 
clouds (with masses < 10"* Mq for R < lOpc) would be larger 



Fueling Seyfert Galaxies 



23 



10 



^ 7 
o 




20 



80 100 

a [km s'^] 



300 400 



Fig. 7. — MBH~cr relation measured by'Tremaine et al." r2002^ (black solid line; Mbh ct''), with estimated intrinsic 0.27 dex scatter (green shaded range), and 
the predicted relation at low-cr from our model for Seyfert triggering as calculated in § 15.11 given an energy-based feedback coupling (dashed line; Mbh ct') 
or momentum-based coupling (dotted line; Mbh «; cr^). Note that these relations apply only at low a where they lie above the solid line. Yellow range shows the 
additional scatter expected (for the energy-based coupling, but similar for the momentum-based coupling) calculated from the size spectrum of molecular clouds 
in § 15.21 For comparison, observations are shown from the compilations of Tremaine et al. 1 2002) (black open circles), Onken et al. 1 2004 ); Nelson et al. ( 2004; 
[Peterson et al. 1 2004) (blue filled circle s), Ba rth et al. 1 2005) (red squares), Greene & Ho ( 2005a) (cyan diamonds), and of POX 52 from Barth et al. 1 2003) 
(purple triangle) and of NGC4395 from lPeterson e t al. (2005) (black star), for which only an upper limit to a exists iFilippenko & Ho 2003). The observations 
tentatively favor the shallower predicted slope and larger scatter at Iow-Mbh from our model, but this remains to be tested with larger samples and more accurate 
mass determinations. 



than the Mbh expected from the co-formation of spheroid 
and bulge only for very small systems a < 25kms~', below 
the range of interest (although these terms do contribute a 
Coul omb logarithm term to the rate of black hole activation, 
see § 14. 1> . However, since the typical timescale for an en- 
counter with a sufficiently large cloud is fevent ~ lO'^yr, up 
to exp(-fH/fevent) ~ 21% of systems may not have had an en- 
counter with a molecular cloud of substantial size at any point, 
and will have an Mbh consistent with the Mbh - cr relation for 
black hole-spheroid formation (Mbh « cr^), making this the 
approximate -Icr lower limit in the Mbh - c relation. 

In Figure we show, in addition to the predicted Mbh-ct 
relation from molecular cloud accretion, the predicted scat- 
ter at each mass (combined yellow and green shaded ar- 
eas). We assume a constant 0.27 dex scatter for the spheroid- 
black hole co-formation Mbh - c mechanism with Mbh oc ct^ 
(|jremaine et al. 2002), shown as the shaded green range. 
We could equivalently show the scatter as a function of 
mass from simulations of merger-driven spheroid formation 
(Robertson et al. 2005), which span a wide range in final 
bulge/spheroid masses and have approximately constant ~ 
0.3 dex scatter 

The observations indicate that the scatter increases at small 
black hole mass. Greene & Ho (2005a) estimate a scatter of 
^0.4 dex primarily in the range cr~ 100-200 kms"', and this 



appears to increase marginally to ~ 0.6 -0.8 dex if we con- 
sider the objects below a ~ 75 kms"'. The increased scatter 
our model predicts at low black hole masses is evident. Quan- 
titatively, however, this change in scatter in the observations 
is only marginally (~ la) significant if we consider points at 
<T < 200km s"', where the predicted (yellow) scatter begins to 
increase over the constant (green) scatter, owing to both the 
relatively small number of points at the lowest masses and the 
relatively large measurement errors. 

5.3. Accretion of Molecular Clouds at Low Mbh 

The condition for blowout, assuming the black hole energy 
or momentum couples to the surrounding medium, implies 
a shallower Mbh - cr relation, with increased scatter, at low 
black hole masses. However, a precondition for this is that 
the black hole be able to efficiently accrete a molecular cloud. 
If the cloud is moving at a velocity cf then the crossing 
timescale during which the cloud interacts with the black hole 
is 

._=-^ = 1.0xl0V(y^)(y^) , (98) 

significantly shorter than the Salpeter time, ts = 4.2 x 10^ yr, 
so the black hole cannot grow by more than a factor of ^ 3 
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in mass in this time interval. Any remaining accretion must 
come from mass captured in the passage of the cloud. 

For a cloud moving at a bulk velocity cf^^, the portion of 
the cloud which passes within a radius 



blowout" region 
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can be "captured" (i.e. bound to the black hole). On average, 
a fraction ~ (R/Ri.{f (given by the cross section of a cloud 
passing over a black hole which does not successfully capture 
it) of the cloud volume/mass is then available for accretion. 
For large black holes, this is more than sufficient to power 
the blowout, and indeed much will not be accreted once the 
blowout begins. For small enough black holes, however, rela- 
tively little mass can be captured. The captured mass is 
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(with a maximum at M^apt = M^) and so a low-mass black hole 
can grow only by at most about one order of magnitude from 
an individual interaction with a typical cloud. If the black hole 
grows rapidly from the beginning of the interaction, this will 
still only increase this by a further factor ^2-3, based on the 
comparison of the crossing time to the Salpeter time above. 
For the black hole masses plotted in Figure0 this is sufficient 
to account for the movement of black holes from the spheroid- 
black hole co-formation Mbh - c relation (Mbh oc a'^) to the 
molecular cloud-dominated Mbh -c relation (Mbh oc ct^), and 
to provide the mass to power the "blowout". Of course, a sig- 
nificant (+lcr) fraction of black holes can interact with larger 
clouds or torii /disk inflows of sizes up to ^ 400 pc (as calcu- 
lated in § 15. 21 above), allowing a factor ~ 200 growth in the 
black hole mass. 

However, for low initial black hole masses '-^ 10^- 1O^M0, 
as for seed black holes from Population III stars, this is only 
marginally sufficient to grow black holes to substantial masses 
(to subsequently act as seeds for further Seyfert accretion or 
growth in mergers). These seeds likely formed at high red- 
shifts z ~ 10-20, when gas densities were large and gas frac- 
tions were of order unity, implying a higher effective "colli- 
sion rate" with dense gas. The Eddington ratio in these phases 
for undermassive black holes is m w 1, so the timescale for 
this growth from the captured mass is « 1 - 2 x 10^ yr, short 
enough compared to the age of the Universe for the buildup 
of seeds which could then grow larger in later mergers. 

5.4. Mass Gain in the Blowout Phase 

Black holes can also grow during the "blowout" phase. For 
a power-law decay of the accretion rate with in = it/t„;)~^'^ 
(Equation 1521 ). the mass gained can be estimated simply if 
it is small and we can approximate M = mMgda = 'wMbh/Js, 
where Mbh is the initial black hole mass. This yields 
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Fig. 8. — Illustration of the geometry at first interaction of a cloud and 
black hole, with cloud radius Scl. black hole radius of influence Rbh, dust 
sublimation radius Sjusi, and cloud velocity Vd ~ c^'*. 



Because liimm <C 1, if jyi > 1 then the above becomes 
AMbh f„-, llL 
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~ tm/ts <C 1, and indeed the mass gained in the blowout is 
small. Similarly, r/£ = 1 contributes an additional logarithmic 
term ln(l/mmin), which still leaves AMbh/Mbh ^ 1- 

For rjL < I, this becomes slightly more complicated, but ul- 
timately gives a similar result. Because mmin ^ 1, this implies 



AMbh/Mbh oc m, 



Equation 



Using the determination of mmin from 



this simplifies to 
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where firji) is a function of rji which varies smoothly from 
e.g. f(T]i = 1/2) = 0.4 to f(T]L ^ 0) ^ 1.8. Therefore, for 
sufficiently low values of rji, the accretion rate decays slowly 
enough that some mass gain is possible in the blowout. 

However, the growth is still relatively small, a factor < 2, 
and less for lower-mass disks. A full numerical calculation 
allowing for the increase in mass with Iow-t^l changing the 
Eddington luminosity as the accretion rate decays also gives 
a small correction of a factor < 2 to the duty cycle at the low- 
est luminosities. This is consistent with detailed simulations 
of X-ray driven wind s, which sugges t efficient wind driving 
for L/Lem > 0.08 (Balsara & KroHk 1993), implying (given 
the large observed energetics of AGN outflows) that most of 
the energy of the system should be generated in this early 
phase. This is also a potential source of increased scatter in 
the Mbh - cr relation, but is at most comparable to th e observed 
scatter (e.g., Tremaine et al. 2002; Novak et al. 2005) and that 
expected from merger-driven co-formation of spheroids and 
black holes, and is much smaller than the scatter expected 
from the accretion of different mass molecular clouds given 
the blowout conditions as discussed in § 15.21 

6. OBSCURATION AND THE CLASSIC MOLECULAR TORUS 

6.1. Obscuration as a Function of Luminosity 

If a molecular cloud is accreted and not all is subsequently 
disrupted by feedback, then some fraction of sightlines to the 
black hole will pass through the surviving portion of the cloud 
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and will be highly obscured. From our momentum coupling 
considerations of blowout, we can estimate this fraction as a 
function of luminosity or black hole mass. 

Consider a cloud approaching the black hole at the time 
where the black hole is about ~ R^i from cloud center (i.e. first 
"touching" the cloud), as illustrated in Figurels] Since /?bh ^ 
Rc\ this indeed defines their "initial" interaction. The cloud 
is moving with a local velocity cf^'', at an angle ip relative to 
the axis connecting the black hole and cloud center Define 6 
as the angle from this axis to a given point in the cloud, with 
the black hole at the origin (defining a ray r), and let <j) be the 
angle about this axis of symmetry. Thus the cloud velocity 
relative to the black hole is v = cf^^v with v = sinip x- cos ipz- 
The linear momentum of a given parcel of gas along such a 
ray originating at the black hole, against the direction of the 
ray, is then 

d|p| = dMcf'^(-r-v) 
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(104) 
and r 



= 2Mci cf'^ cos ^ cos'* dcos , 
where the third equality comes from integrating over 
(from r ^ Rbh ^ ^ci to r = 2Ra cos 9). 

Consider the portion of the cloud with > 6*0, then the total 
momentum input to "clear" a sightline in this fraction of the 
cloud is 

\p\ = ^ cosV'Mcicf'^cos^eio. (105) 



Since by our definitions < -0 < 7r/2, {costp) =2/tt, so 

cos^6'o. (106) 



4.77- 



If the radiation from the black hole is directed isotropically 
and impacts the cloud, at least initially, via momentum cou- 
pling, the momentum imparted in some time interval At is 
then dp = 77(L Af/c)(dfi/47r), where dft is the opening an- 
gle considered and 77 is the coupling efficiency. Considering 
the range above, with Oo < < tt /2, this gives Aft = J dft = 
In cos 9o, so the total momentum input from the black hole in 
this region is 

1 LAt 

P^n = ^V cos 6*0. (107) 

2 c 

This will be sufficient to overwhelm the momentum of the 
cloud and blow out material with > 0q for 
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where we have used Af = t,,-, ^ Ra/o', the characteristic 
timescale for the blast wave impact. 

Thus, for sufficiently low luminosities or black hole masses, 
there is some solid angle in which feedback cannot "punch 
through" the cloud, and the cloud will remain as an obscurer. 
The covering angle for this is 

/=i(l-cos0o). (109) 




Fig. 9. — Predicted Seyfert 2 fraction as a function of luminosity, calculated 
from our modeling of a pressure-driven blast wave (bl ack li nes). Dashed line 
uses the simplified luminosity criterion of Equation 1108L dot-dashed line 
includes the more detailed dependence on Mbh (integrating over the high- 
accretion rate period), solid line further includes integration over a spectrum 
of cloud sizes. Points show the observed Type 2 fraction from the local AGN 
luminosity functions of iHao et ali (i20Q 5). using Schechter (circles) and dou- 
ble power-law (squares) function fits to Ha (blue), [Oil] (red), and [OIII] 
(green) observations, converted as in Figure |2] For comparison, the fitted 
Type 2 fraction from Simpson |2005), determined from the samples therein 
and from Ueda et al. |2()03); Hasinser 1 2004,) :. Grimes, Rawlinss, & Willotil 
( 2004), and dominated by quasars with < z < 3 (i.e. merger-driven quasar 
activity and obscuration in our modeling; see Figure |6| and Hopkins et al. 
2006a) is shown (cyan line). 



This is similar to our deriv ation of the blowout criterion as a 
function of mass in § 15.11 and here we incorporate the geo- 
metric effects of "partial" blowout. These partially obscured 
objects will be visible as Seyfert 2's, with properties consis- 
tent with the canonical molecular torus unification model. 

In Figure 121 we show the predicted fraction of Seyfert 2's 
as a function of l umino sity from the luminosity-dependent 
form of Equation dlOSt . For comparison, we plot the esti- 
mated Seyfert 2 fraction from observations as indicated, in 
the manner of Figure|2] The predicted scaling agrees with the 
observed trend, and also defines a cutoff at L ~ 10" - lO'^L©, 
i.e. Mb ^ -19 to -21, above which the requirement for obscu- 
ration is cos 6*0 > 1; i.e. the black hole is sufficiently luminous 
to destroy the entire cloud regardless of angle. 

There are, however, two detailed effects which will smooth 
this cutoff. First, the pure luminosity form of Equation jl08l l. 
while giving the appropriate scaling, is not strictly appropriate 
since the blowout in a given direction may not depend on the 
instantaneous Seyfert luminosity, which may have decayed to 
quite low luminosities at later times, but on the momentum 
injected in the initial, high-luminosity "blo wout" phase. This 
is given by re-considering Equation dlOSt . as a function of 
black hole mass, with m = 1 for a period f,„ as determined in 
§ 13.21 (technicallv we integrate from ni = 1 to m < 0.1, since 
we define a differential time per logarithmic interval in ac- 
cretion rate). This then provides a more precise estimate of 
the fraction which will have been "blown out" by the blast 
wave. Accounting for this yields the dot-dashed line in Fig- 
ure |9] At a given luminosity, there is a distribution of black 
hole masses at various accretion rates, and this effect smooths 
the obscured fraction as a function of luminosity. However, 
this still defines a characteristic Mbh ~ lO^Msim ^ Md above 
which the blowout will completely disrupt the cloud, with a 
corresponding Eddington luminosity L ~ lO'^L©. 

The second effect is the spectrum of cloud sizes. A more or 
less massive cloud also implies a correspondingly smaller or 
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larger covering angle of material which cann ot be expelled. 
We use the spectrum of cloud sizes in § 15.21 to calculate the 
dispersion in the Mbh-o" relation expected from this mode of 
accretion, and an identical calculation gives the dispersion in 
Mci and correspondingly, via Equation ( I108> . the variation in 
obscured fractions as a function of Mbh- Averaging over this 
distribution yields the total obscured fraction as a function of 
Mbh and L, shown as the solid black line in Figure|9] 

These more detailed effects produce significant differences 
at intermediate masses, where the transition between ob- 
scured and unobscured populations is smoothed out. The 
agreement with the observations is improved, suggesting that 
this momentum balance argument provides a plausible esti- 
mate of the obscured fraction as a function of luminosity. 

Note also that our argument does not necessarily depend on 
the black hole luminosity affecting directly the surrounding 
medium via radiation pressure or other momentum couplings. 
Rather, because the blast wave impacting the medium will be 
pressure-driven, 77 as used above represents an effective effi- 
ciency of the Seyfert energy and momentum "ending up" in 
the initially driven blast wave by the time it impacts on the 
surrounding medium. 

Our prediction is qualitatively similar to those of various 
modified "luminosity-dependent" torus models, and the fun- 
damental qualitative point, that black hole feedback which is 
stronger in higher-luminosity systems is able to clear a larger 
opening angle, still obtains. However, our result is distinct 
from some of these models in several respect s. For example, 
the "receding torus" model (e.g., Lawrencel ll991h assumes 
that the opening angle increases as the inner torus radius in- 
creases, owing to a larger dust sublimation radius rdust L'/^ 
with increasing dissociating luminosity. However, while this 
gives a reasonable scaling of the obscured fraction for a geo- 
metrically flat (i.e. rectangular torus cross-section) torus, it is 
clear from Figure|8lthat for an increasing torus radius with r 
(i.e. a circular torus cross-section or elliptical/spherical cloud, 
and conical opening angle models), removing a spherical re- 
gion around the origin at the black hole has almost no effect 
on the covering angle of the cloud. This is true until rdust ^ ^ci 
(for example, even at rdust = ^ci, for a spherical cloud or circu- 
lar torus cross section, this only decreases the Type 2 fraction 
by a factor of two from rdust = 0). In other words, the dust 
sublimation radius would have to reach ~ 10- lOOpc before 
this mechanism would explain a scaling of obscured fraction 
with luminosity in this fueling scenario, into (or even beyond) 
the traditional narrow line regions. 

Qualitatively, our mechanism is more s imilar to th e 
disk-w ind model of Koni gl & Kartil (11994 : lElvisI J2000I) : 
iProga & Kallman (2004), in which radiation pressure flattens 
the initially vertical outflow and dust distribution in objects 
with higher bolometric luminosities. However, it can oper- 
ate on much larger scales and depends on momentum balance 
from the blast wave (not necessarily direct radiation pres- 
sure) ejecting material within some opening angle, rather than 
merely flattening the distribution. It is also worth noting the 
striking si milari ty of the obscured fraction calculated from 
Equations ( I108> and ( I109> and the functional form of empir- 
ical calculations which allow for e.g. evolution of the inner 
torus radius and torus height as a function of luminosity (e.g., 
ISimpsonll2005h . suggesting that feedback balance does cap- 
ture the critical dependence on black hole and fueling proper- 
ties. 

For comparison, the fitted Type 2 fraction as a function of 



luminosity in higher-redshift quasars fro m ISimpsonI 1 20051), 
determined frorn the SPSS and from lUedaetaLl (20031) : 



iHasinged ( 12004 : iGrimes. Rawlings. & WillottI ( 12004 7 is 
shown as the cyan Une in Figure|9] While a similar qualitative 
trend is apparent, the Type 2 fraction is much higher, com- 
pletely dominating at low-luminosities and contributing a sub- 
stantial ^ 30% portion of the population even at Mb < -23, 
whereas in the local samples the Type 2 fraction is at most 
comparable to the Type 1 fraction at low luminosities and neg- 
ligible (^ 10%) at high luminosities. 

The cause of this difference is suggested by Figure |6l in 
which the Seyfert (i.e. quiescent, cold gas driven) and quasar 
(i.e. merger-driven) luminosity functions are each plotted at 
several redshifts. Even by z ~ 0.1-0.2, the quasar popula- 
tion is larger at these luminosities, and indeed the Simpson] 
(2005) sample is drawn from surveys which are dominated 
by quasars at high redshifts spanning the wide range 0.15 < 
z < 3. These objects, in our modeling, are triggered by differ- 
ent processes, and therefore should not necessarily follow the 
same trend - indeed, the discrepancy between the two obser- 
vations cautions against extrapolation of obscuration models 
(e.g. geometrica l models) cahbrated for local, non -interacting 
populations. In iHopkins e t al.' (>2005a^b^. '20 06a^ we discuss 
this obscuration in a merger-driven context in detail, an d 
show that the trends of "Ueda et al. (2003); Hasinger (2004"); 
Grimes. Rawlings. & Willott (2004); Simpson (2005) can be 
reproduced when the primary source of obscuration is pro- 
duced by larger scale (> 50 pc) gas inflows from the mergers 
which themselves power accretion, and the distribution of col- 
umn densities derives primarily from time-dependent rather 
than line-of-sight dependent differences. A similar qualita- 
tive trend is found because bright systems are more likely to 
be near their peak luminosity, when they are able to expel ob- 
scuring gas and dust, which is similar in principle to the model 
we have described here for obscuration of cloud-fueled AGN. 
However, in our model for non-interacting Seyferts, the ob- 
scuration is removed along certain sightlines at the beginning, 
not at the end of black hole growth, and the distribution of Nh 
does, presumably, derive from geometric effects, as in canon- 
ical torus or disk-wind models. 

6.2. Building a Torus 

Although it is not a requirement of this model, our picture 
for AGN fueling can both derive from and seed the tradi- 
tional molecular torus of standard AGN unification models. 
Assuming that a significant fraction of an accreted cloud is 
captured on passage near the black hole, essentially requiring 
Rci < GMbh (cf'^T^, or Mbh > 10^ Mq (ensured for the black 
hole masses of interest and for the predicted Mbh - f relation 
for this fueling mechanism, see § |5|i, then it will presumably 
develop an orbit, circularize, and some fraction of the cloud 
will be tidally disrupted. 

Given an initial impact parameter b and velocity cf^^, the 
material can settle into a Keplerian orbit determined by con- 
servation of angular momentum. For the bulk of the cloud, 
b ^ Rci^ Rbh, so the circular velocity is determined by the 
spheroid potential. 
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where a is the bulge scal e-len gth. Using a = 
10kpc(cr/200kms"')^ from § 13.11 with the appropriate 
factor of <j) to convert from observed cr to a physical Vc 
and a assuming a Hernquist (1990) spheroid profile yields 
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Vc w 200 km s ' \J r/kpc. Since this will include most of the 
mass, demanding the final angular momentum ^Mrvc (in 
detail = (3/1^/4+ X)MrVc for a torus with h = h/r) conserve 
the initial ^ Mj^ cic?" *^ (again, in detail M R^\cf^ simj; for 
the geometry in § 16. U yields a torus radius 
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independent of a, and assuming h ^ I as we show be- 
low. Thus the characteristic torus size < lO's of pc (up 
to a maximum ^ lOOpc) inferred from a wide variety of 
observ ations (Krolik & Beaelman 1988; Granato et al. 199^ 
Bock et al. 2000; Schinnerer et al. 2000; Galliano et al. 2003; 
Radomski et al. 2003; Weiaeh et al. 2004; Ja ffe et al. 2004; 
PrietoetalJl200 4; Elitzur 2005) and simulations of torus 



structure and stability with radiative transfer models (e.g., 
|Pier & Krolik 1993; Granato & Danese 1994) is a natural 
consequence of the scales of molecular clouds, and is "built 
into" this picture of Seyfert fueling. 

Given this radius, whether the molecular cloud is circular- 
ized into a strict torus with volume 2tt^R^ (assuming h ^ I) 
or remains spherical/ellipsoidal, the change in e.g. density and 
other gas properties is well-defined (varying by e.g. 20% us- 
ing the exact solutions above and assuming it goes from a 
sphere to a torus). Therefore, the characteristic density is 
also given by that of the original c loud, enhanced by the com- 
pression expected from Equation il 1 II above, explaining the 
molecular and dust densities and masses inferred for molecu- 
lar torii around Seyferts, with « ~ 10^- lO'^cm"-'. This gives 
a Une-of-sight column density through the torus plane of 

A^Hi = 2/?,«t«2/?cinci (^)^» 1.6 X 10^^ cm'^ (112) 

This is similar to the column densities of torii inferred 
from phenomenological fits to the com plete hard X-ray col- 
umn density distributions of Seyferts ( Urrv & Padovani'1995'; 
[Treister et al. 2004; Mainieri et al. 2005) and from synthe- 
sis models fitted to the X-ray backgroun d (Setti & Woltie: 
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[1^89; Ma dau et alJ fT994l IComastri et alJ ri995,: .GilU 
Il999, 2001i lTreister & Urryl2005l) 

We can also calculate the expected height or flattening of 
this torus, as a function of radius and Seyfert properties. If 
the molecular cloud, with density n^i ~ lOOcm""* collapses 
from the ISM with mean density hism ^ 1 cm"-', and magnetic 
flux (= BdA ^ BR^) is conserved, the magnetic field will be 

enhanced by Bci/^ism ^ ("ci/niSM)^''''- For typical magnetic 
fields providing pressure support in molecular clouds 
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with Bel = BiSM (nci/niSM)^/^ and Md = (47r/3)/?^, p. For a dif- 
fuse ISM BisM ~ a few /iG, this yields 



(114) 



giving the standard self-consistent picture for the magnetic 
support (in addition to support from turbulent motions) for 
molecular clouds, with these field strengths characteristic of 
the disordered fields of the Milky Way ISM (e.g. Troland & 



Heiles 1986; Rand & KuUcarni 1989; for a review see e.g. 
Ferriere 2001). Increasing magnetic field strengths are also 
observed towards the center of the Galaxy, up to a factor of 
^10 larger, supporting higher densities and even larger giant 
molecular cloud complexes with Md up to ~ lO^Mp, (e.g., 
[ Rand & Lv ne 1994; Mvers et al...l995: Mvers & Khersonskvl 
ll995UCrutch er 1999). 

Because capture and circularization does not completely 
erase the properties of the molecular cloud, the magnetic 
fields on cloud scales can be preserved, continuing to pro- 
vide pressure support. This pressure defines an effective 
sound speed c| = (B^[/47r)p"', which for the above is cb ~ 
20kms"' (BiSM/4/iG). Note that any other source of a sim- 
ilar pressure which is non-thermal (i.e. does not scale adia- 
batically), such as e.g. turbulent motions, will have a similar 
effect. Since this dominates the pressure support of the torus, 
and using the circular velocity at Rt ^ Rbh (i e- dominated by 
the spheroid potential) the scale height is 
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with a maximum of Ht ^ Ra- Thus, the torus is not ex- 
pected to be strongly compressed, at the outer radii, with 
Ht ^ Rt, again similar to traditional phenomenological mod- 
els and observat ional estimates of the torus structure and 
covering angle ( Krolik & Begelman 1988; Antonucci 1993^ 
Urry & Padovani 1995; Risaliti et al. 1999; Schmitt 2001). 

The timescale for all of this to occur will be of order the or- 
bital timescale, foit ^ f/'^c, or ~ Ra/cf^^ ^ 10^ yr for the bulk 
of the cloud at Ri ^ R^i. Thus, if the cloud is captured, it will 
begin to reach an equilibrium in the blast wave and orbit with 
the circular velocity in roughly a crossing time. It is there- 
fore expected that it will be seen as such over the majority of 
the Seyfert lifetime. There is a period of time ^ 10^ yr, com- 
parable to the timescale at highest accretion rates m > 0.1, 
during which the system may be settling into equilibrium. We 
have shown the transition from cloud to torus implies little 
change in the cloud properties, so observational signatures of 
this transition phase, if it occurs, may not be obvious. 

7. GLOBAL CONSEQUENCES FOR THE HOST GALAXY 

Given the effects of the black hole feedback we have mod- 
eled, and the fact that the blast wave can propagate to large 
(~kpc) scales on reasonable (~Gyr) timescales, in principle 
there could be a significant impact on the host galaxy ISM or 
gas structure. The black hole feedback expected from growth 
in mergers can sweep a large fraction of the remaining gas in 
the galaxy into an unbound blast wave, heating the remaining 
jas to the virial temperature and terminating star formation 



gas to 



llHopkins et al. 2006a c; S pringel et al. 2005a). However, the 
case of fueling by molecular cloud accretion is quite different, 
and we find that the blast wave is much weaker than required 
to significantly impact the host galaxy. 

If the blowout is defined by the criterion that Md be accel- 
erated to v'esc ^ (i.e. that the local gas be unbound), then 
because Md <C Mism, a collision between the blast wave and 
the ISM will at most accelerate the ISM to ^ (Md/MisM)^, 
if momentum considerations are most important. For typi- 
cal Galactic scales this gives a bulk speed <C Ikms"', which 
is completely negligible on the scales of the ISM. Even 
under the extreme assumption that all the blast wave en- 
ergy can be converted into bulk motion gives a velocity ^ 
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\/Mc\/MisM <7 ^ Ikms"'. Thus, the blast wave itself is not 
expected to significantly impact the host galaxy ISM. 

If, instead of the momentum or energy in the blast wave 
shell, we make an even more extreme assumption that all 
of the feedback energy or momentum flux from black hole 
growth impacts the galaxy ISM, the energy input could be 
much larger We further make the maximal assumptions that 
rji is smal l, so the black hole grows substantially (by a factor 
~ 2, see § I5.4> . and that the affected ISM is unable to radiate. 
Let us also consider the case of direct energy coupling (i.e. 
conversion of all feedback energy into bulk motion), which 
gives an increased impact of feedback by a factor ^ cja ^ lO-' 
over e.g. momentum-based coupling. Then, the feedback en- 
ergy from the black hole is ~ 77 e, MBH Ai7, where lS£t is the 
covering angle of the disk to the black hole, since the radia- 
tion is isotropic. Note that collimation will tend to be perpen- 
dicular to the disk and along the angular momentum vector, 
meaning it will only moderate the impact on the host galaxy. 

The binding energy of the disk gas is ^ /gas^^dt^'?- Using 
the fact that Mbh ^ where Mb is the bulge mass, this 
implies that the black hole feedback will be significant (i.e. 
comparable to the binding energy of the disk ISM gas) when 
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/gas Mi 

Since the disk covering angle (weighted over the mass/volume 
of the disk) is 9 ^ hd/r^ ^ cj/vc, and the covering factor is 
(1 -cos6l) ^ 6*72, and we expect ^ - 0.001, rj - 0.01, cf'^ - 
lOkms"', we can solve for the maximum Vc for which black 
hole feedback will have a large impact. This gives 

' f \ -1 

./gr- 
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^o.ooiy V 0.1 yj 

i.e. Md < 10^ M0 with Mbh < 10^ Mq. Again, this assumes 
the entire energy output of the black hole has a chance to cou- 
ple to the ISM with no allowance for that ISM to dynami- 
cally relax or cool, and thus represents the most extreme es- 
timate of which galaxies could be affected. Such low-mass 
galaxies will typically be later types than the Sa/b's with typ- 
ical B/T^ 0.1 above, more likely < 0.02 for Sc/d and Sm/Im 
galaxies, and gas fractions potentially quite a bit higher than 
^0.1, making the host even more robust against black hole 
feedback. If we consider a momentum coupling (^ rjerMgiiC 
compared to /gasMj v^) criteria instead, then we obtain 

.o.i/vo.i/ vi.o/ 

\ /B/TmI/3 
.0.001/ V 0.1 
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a negligible effect. 

8. THE CONTRIBUTION OF STELLAR WE^D FUELING 

For comparison, we briefly consider the contribution of 
stellar winds and hot (virialized) gas fueling to AGN activ- 
ity. First, consider the contribution of stellar winds from the 
galaxy as a whole - these will, in general, shock against one 
another and virialize, contributing to the hot gas reservoir 
which can be accreted through the Bondi radius. Assuming 
the gas is virialized gives = 3(f) a^, where the dimensionless 
factor (j) depends on the density profile and converts between 



the observed mean projected velocity dispersion within the 
effective radius (a) and the potential at z = ((/> w 10.1 for a 
Hernquist [1990] spheroid profile). We further use the radial 
hydrostatic equilibrium condition to solve for p, which gives 
a profile having a core with a characteristic njsM ~ 1 cm""* at 
the center Using this and the Mbh ~ c relation, we obtain 



m = Ana 



G^MBHtsP ( o 
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where we have also taken 7 = 5/3 (giving a = 1 /4 for the 
B ondi-Hoyle solution) in the second equality. 

Ciot ti et all lfT991l) consider a detailed numerical calcula- 
tion of injection of gas into the hot ISM from steady stel- 
lar mass loss, supernova events, and galaxy-scale inflows and 
outflows. For each of their range of conditions, we consider 
the temperature, gas velocity, and density at the trans-sonic 
radius for inflow (in the solutions they consider with a central 
mass concentration) and compare the implied Bondi rate with 
that estimated above, and arrive at a similar result (to within a 
factor ^ 2-3). This also agrees with the Bondi rates estimated 
from systematic measurements of hot gas temperature and 
density profiles in the inner regions of a number of quiescent 
ellipticals (e.g., Pellegrini 2005; Soria et al. 2005b). The char- 
acteristic rh ~ 10"^- lO"** from galaxy scale hot virialized gas 
appears to be a robust result. Note that while the total stellar 
mass loss rate of a quiescent elliptical can be ^ 1 Mm yr~' un- 
der typical conditions, the calculations of lCiotti etaUiMl 
and others have shown that a far smaller fraction of this gas 
actually reaches the center of the galaxy, and it does so either 
at super-cr bulk velocities (if un-shocked) or ~ a sound speeds 
(if shocked and virialized), giving a much lower actual Bondi 
rate for black hole accretion. 

Following Soria etal. (2005a b), we also estimate the to- 
tal stellar mass loss within the black hole radius of influence. 
Winds within this radius will be captured by the black hole, 
and be available for accretion. Note that, strictly speaking, 

this is appropriate only for wind velocities < \[^<^ (the 
spheroid escape velocity at the center), as otherwise the cap- 
ture radius is smaller than the radius of influence. For typi- 
cal slow wind speeds ^ 100- 300kms"', however, this is at 

most comparable to typical ellipticals \J^a, making this a 
relatively small source of uncertainty. This will, however, re- 
strict the ability of e.g. a young stellar cluster near the nucleus 
(with much higher M in winds than older stellar populations) 
to contribute significantly to moderate or high-;7t populations, 
as with typical v„, up to ^ 1000 kms"', the effective volume in 
which such systems can have their winds captured is dramati- 
cally reduced, and furthermore the large bulk velocity reduces 
the estimated Bondi rate. 

For steady mass loss from o lder stellar populations, we 
adopt the estimate of lCiotti et a l. ( 1991 ) from population syn- 
thesis for the mass loss rates as a function of the spheroid age 



Gyr. 
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Using a constant mass-to-light ratio A = 
(Mq/Lq g)/(Mbiiige/iB), and modeling the spheroid stars 
with a Hernquist (1990) profile to calculate the mass fraction 
within /?BH = GMbh/ct^, then using the black hole-bulge 



Fueling Seyfert Galaxies 



29 



mass correlation of |m arconi &Hunll ll2()0l to cancel the 
dependence on MBn/Mbuige, we obtain 

'MfiH/Mbulge^ 



1.1 X 



10"^ A I 



-1.3 



(121) 

0.001 yVlOGyr/ 
which, for reasonable values of A (generally below unity since 
these stellar populat ions are old) from the stellar population 
synthesis models of iBruzual & CharlotI ( 1200 3 ) gives m in the 
range ^ 10"^ (for relatively young ellipticals with ages ^ 2- 
4Gyr) to ^ 10"^ (for older ellipticals with ages ^ lOGyr). 

In either case, fueling from stellar winds and quiescent 
Bondi accretion of hot gas are roughly comparable, within 
the range ~ 10"^ - 10"^, without a strong dependence on host 
galaxy properties. Given our modeling of quasar formation in 
mergers with a subsequent decay in the quasar light curve de- 
termined by the evolu tion of the feedback-driven blast wave 
jHopkins et alJ 12006c). we then expect that ellipticals will, 
shortly after forming, be at moderate, declining accretion 
rates for some time, contributing to the -18 > Mb > -23 end 
of the predicted elliptical Seyfert population (i.e. the bright 
end of the local Seyfert luminosity function, but faint end 
of the merger-driven quasar luminosity function). Ho wever, 
on ce a s ufficiently low accretion rate given by Equations (^^J 
& ( I121> is reached, the hot gas in virialized equilibrium and 
stellar mass loss of aging stellar populations near the black 
hole sustains a nearly constant low-level accretion rate. 

Accretion rates of black holes in early-type galaxies will 
then decay from the redshift when the elliptical formed to 
these Eddington ratios and "pile up" at m ^ IQ~^. This 
is supported by direct estimates of the Eddingto n ratio dis- 
tribution in quiescent "dead" elUpticals (^ e.g., iHol l2002t 
iMarchesini et al .''2004': Hopkins et al. 2005e). Furthermore, 
lower mass ellipticals which form later (e.g. Cowie et al. 
1996; see Hopkins et al. 2006b for details of our calculation) 
may not decay quite to these lowest m values, ultimately over- 
lapping in luminosity. Therefore, a large fraction of the ellip- 
tical population fueled by such mechanisms will be found in 
a relatively narrow range in accretion rate and luminosity. 

To estimate the contribution to the luminosity function, we 
must account for changes in the radiative efficiency (e.g., 
iNaravan&Yi 1995; Esin. McClintock. & Narayan 1991 
iQuataert 2001; Naravan 2004) as a function of luminosity and 
the fact that such low-«i systems are observed to be ac creting 
substa nti ally below the Bondi rate ( Fabian & Canizgrgj 
1988| lElandford & Begel manI IT999; .Pi Matteo et al 
200a iNar avan et al. 2000;' lOuataert & Gruzinovl 120001: 
Di Matteo et al. 2001; Loewenstein et al."200r: 'Bow er et'all 
2003; Pellegrini 2005). 
efficiency 

r 0.1 



We do so by adopting a radiative 



0. 



1 4^- 

merit 



if ffi > merit 
if m < merit 



(122) 



a transition between radiatively efficient and ineffi- 
accretion at nic-it ~ 0.01 as suggested by obser- 



for 

cient accretion at merit 

vations of black hole binaries ('Maccarone' '200 31) and 
AGN Eddin gton ratio distributions (Marchesini et al.' 1 2004 
Ijeste3l2005h. and theoretical extensions of accretion mod- 
els (e.g., iNarayan & Yi 1995; Esin. McClintock. & N aravaa 
[T597; Meyer. Liu. & Meyer-Hofmeister 2000). This choice 
for the efficiency follows from ADAF models (Naravan & Yi| 
[l995), also allowing for large m ass loss through winds giv- 
ing sub-Bondi accretion rates (Blandford& Begelman 1999; 
ISoria et al. 200 5b) and broadly accounts for observations of 
local quiescent objects. 



For a characteristic m 10 (i.e. L/LEdd ^ 10 ^) and early 
type black hole masses ^ 10^- IO'^Mq with redshifts of for- 
mation z ~ 2, this implies a narrow, steep distribution in Mg 
peaked around -11 < Mg < -13.5. In detail, we again use 
the determination of the elliptical m distribution and lumi- 
nosity function in Figure |2] from our modeling of black hole 
and spheroid co-formation in Hopkins et al. ( 2006a b) and 
the decays of those light curves in our blast wave solution 
from Hopk ins et al. (2006c), but impose the larger of Equa- 
tions dl 19l l & ( I121> as a lower limit to the accretion rate (for 
f > 1 Gyr) decay. Convolved over the early-type population, 
this gives the dot-dashed line in Figure |2] Indeed, it is sug- 
gestive that the estimated contribution from this mode of ac- 
cretion becomes important at about the luminosity where our 
prediction for stochastic cold gas accretion flattens and may 
begin to fall below the observations. As discussed in § 14.41 
this more involved calculation ultimately implies a relatively 
basic point, that there are a large number of "dead" quiescent 
early-type galaxies with steady-state m ^ 10"^- lO"'* accre- 
tion rates, which make up a substantial portion of the AGN 
and LINER population at the lowest luminosities. 

9. DISCUSSION 

We have developed a model for the light curve and accre- 
tion rate evolution of accreting systems in feedback-driven, 
Sedov-Taylor type blast waves. Our formalism is applicable 
to a wide variety of systems with different equations of state, 
external gas profiles, physical scales, and feedback coupling 
mechanisms. We have applied this general model to the spe- 
cific case of supermassive black hole accretion via "stochas- 
tic" (i.e. not cosmologically-induced) collisions with molec- 
ular clouds (or the inflow of such clouds from disk, bar, or 
torus processes), in quiescent systems with some supply of 
cold (rotationally supported) gas. Feedback from accretion 
energy rapidly unbinds nearby gas and the system is described 
by our feedback-dominated decaying light curve solution. 

Predictions from this model are consistent with many prop- 
erties of low-luminosity AGN, which are indeed observed to 
be in quiescent systems with cold supplies of gas, and our 
picture has many testable consequences. Because the fuel- 
ing mechanism is not cosmological in nature, these predic- 
tions are essentially a priori, and do not require a detailed 
cosmological modeling of the evolution in galaxy properties. 
Moreover, this provides a context for considering such objects 
and their large-scale fueling mechanisms, and for contrasting 
them and their evolution with that of bright quasars driven by 
the cosmological processes of interactions and galaxy merg- 
ers. 

9.1. Comparison with Observations 

Our model reproduces numerous observations of low- 
luminosity AGN, including: 

• Duty Cycles: The duty cycle at high accretion rate s m > 
0.1 is expected to be ^ 1% (compare e.g., jK auffman n et alJ 
I20033IY.1 et alJl2?MlDong & De Roberti J200,5.). This rises 
to imply that a large fraction of galaxies host a low-luminosity 
AGN, with a high active fraction of ^ 20-50% by Mb >-12 
(depending on the location of the cutoff) (e.g., iHao et all 
[2005: Best et al. 2005). The inflows of cold gas observed in 
iLauer et aL ti2005i) in the central regions of low-luminosity 
AGN also appear to follow a similar duty cycle and period- 
icity to our predictions. Note, however, that technically these 
are theoretical upper limits to the duty cycles, for if accre- 
tion proceeds intermittently (i.e. in short, potentially super- 
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Eddington "bursts") the same average accretion rate on the 
timescales relevant for our calculations is maintained, al- 
though the timescale for such bursts is still constrained by the 
observed episodic quasar lifetime (see e.g.. Martini 2004). 

• Seyfert Luminosity Function: We predict a local Seyfert 
luminosity function which agrees with that observed (e.g. , 
iHao et al Im&i luivestad & Holl2()()!l iHuchra & Burglll992h 
over nearly ten magnitudes, -14 > Mb ^ -23. This re- 
sult depends only weakly on the theoretical uncertainties in 
our blast wave model and more uncertain duty cycle calcu- 
lation. Although our prediction may fall short of the ob- 
served luminosity function at low luminosities Mg > -14, this 
is where contamination from star formation becomes a se- 
rious observational concern (see the discussion i n Hao et al. 
12005'; Kauffmannet al. 2003a: Ke wlev et al.l f2001). and also 
where the contribution from quiescent, relaxed ellipticals 
fueled by hot gas accretion and s tellar winds may begin 
to do minate (see § |8] above and ^Ci otti & Ostriked 119971 
1200 U [Pellegrini 2005; Soria et al. 20051^. At high luminosi- 
ties, our predictions map neatly onto the quasar luminosity 
function observed (e.g., Bovle et al. 2000; Ueda et al. 2003; 
[Richards et al. 2005) and predicted from merger-induced ac- 
tivity ( Hop kins et al. 2006a), with an different fueling mech- 
anism but similar feedback processes regulating the decay of 
the light curve. 

• Morphological Properties: From our analysis, we pre- 
dict the contributions of various morphological types to the 
luminosity function. At high luminosities, there will be 
a significant contribution from relaxing ellipticals decaying 
from recent (~ Gyr) merger and starburst activity as seen in 
e.g. Kauffmann et al. (2003a); Sanchez & Gonzalez-Serran3 
(l2003l) : iSanchez et all (|2004); Vanden Berk et al. (2005). In 
the range -20 > Mb > -23 there will be a comparable con- 
tribution from SO and Sa/b systems at moderate to large 
accretion rates, with a potentially non-negligible but rela- 
tively small contribution from Sc/d systems (and smaller 
bulge/black hole systems). These Sa/b and SO systems dom- 
inate the low-luminosity activity over the range -15 < Mb < 
- 20. Both of these tre n ds are seen observat io nally, e.g. 

fauffmann et all (|2£)03a'); 'Sanchez et al.' ('2004'); "Best et al.' 
005) and specifically by Dong & De Robertis (2005) who 
give a detailed breakdown of AGN activity in different late- 
type morphological types. These objects are not necessarily 
associated with mergers or interactions, unlike the highest- 
luminosity quasars (not seen at z = 0), again consistent with 
the observed morphologies (e.g., Kauffmann et al. 2003a). 
Below these luminosities (i.e. luminosities characteristic of 
the lowest luminosity AGN and LINERs), "dead" early- 
type galaxies fueled by e.g. steady accretion of hot (viri- 
alized) gas and stellar winds from old stellar populations 
make up a l arge fraction of the ac ti ve population, as seen by 
IHo (2002); Heckmanetal. (2004); Marchesini et af. (2004j|); 
Ijester (2005); Pellegrini (2005) ; Soria etal. (2005b). Inde- 
pendent evidence for our hypothesis comes from observa- 
tions of Seyfert clustering, which suggest that Seyfert hosts 
have masses of typical gas-rich, late-type (typically low-mass) 
systems, whereas large black holes with low accretion rates 
show clustering typical of l arger mass, gas-poor ellipticals 
(fc:onstantin & Vogelevl2006ft . 

• Outflow Properties: According to our model, even low- 
luminosity AGN should be associated with some locally 
driven Sedov-Taylor like wind or outflow, although this will 
not necessarily globally impact the host galaxy. Many ob- 
servations find such an association, with characteristic en- 



trained masses ^ 10''- 10^ M0 and small-scale wind dynam- 
ical times ^ 10'*-10^yr, driving outflows with mass out- 
flow rates several times the AGN accretion rate and implied 
energy coupling efficiencies ^ 0.01 - 1 (see Veilleux et al] 
2005, for a review), consistent with feedback-driven de- 
struction of stochastically accreted gas clouds. The velocity 
and t emperature structures (e.g. , Shopbell & Bland-Hawthorra 
T998|lCrenshaw et al.'200tf.'&olik & Ki-iss 2001: Kaspi et all 
20021 iRice et al. 2006, and references therein), small en- 
ergy injection zones (e.g.. Smith & Wilson 2001), and de- 
tection rates (e.g.. .Crenshaw et al. 1999; Rupke et al. 2QQ§) 
of observed AGN wind s, even in obscur ed and starburst- 
dominated systems ("e.g. jRupke et alJ200 5^. support this pic- 
ture. Of course, it has long been recognized that these 
observations can be modeled by Sedov-Taylor blast waves 
driven from near the black hole (e.g., Begelman et al. 1983 
Shlosman, Vitello, & Shaviv 1985; Konigl & Kartie 1994 
.Murrav et al.. .1995.: .Elvis, .2000,) . but our modeling self- 
consistently links the outflow and accretion rate evolution. 

• Eddington Ratio Distributions: We predict the distribu- 
tion of Eddington ratios in low-luminosity AGN, as a function 
of e.g. velocity dispersion, black hole mass , and lu minosity, in 
agreement with that observed bv lYu etal.l(l2?)?)5ll . over three 
orders of magnitude in m, for a range of cr ~ 70- 120 kms"'. 
This includes both the power-law behavior of the Edding- 
ton ratio distribution set by our blowout solution and the 
turnover at the lowest rii determined by the finite time be- 
twee n new "excitati on" events, resolved from comparison 
with lYu et alJ (l200.5h at the 4 - 6 cr level. Predictions in- 
cluding the distribution of post-merger or hot gas-fueled el- 
lipticals also agree well at larger a, but this is less directly 
related to our modeling. The cumulative Eddington ratio dis- 
tribution, combined with that predicted for "dead" ellipticals 
from Hopkins et al. (2006a), m ay be bimoda l, as suggested 
by e.g. Mar chesini et al. (2004); Jester (2005'). However, this 
bimodality is easily washed out by relatively small changes in 
our estimated duty cycles, and is therefore not as strong pre- 
diction of our model as is e.g. the shape of the Eddington ratio 
distribution in active late-type systems alone. 

• Eddington Ratio vs. Luminosity and Mass: We find that 
the mean or median Eddington ratio is a strong function of 
the observed (intrinsic) AGN luminosity. In other words, the 
observed low-luminosity Seyfert luminosity function spans a 
wide range in Eddington ratio and a relatively small range in 
black hole mass, so the observed luminosity function is pri- 
marily an Eddington ratio sequence determined by feedback- 
driven blowout. This a grees with e.g. IHao et alJ (|2005); 
IVanden Berk'etall (l20()5l) . who find little correlation between 
AGN and host galaxy luminosities, suggesting that the ob- 
served luminosity is primarily a function of Eddington ratio 
with characteristic black hole masses ^ 1O^M0. The mean 
Eddington ratio is also a function of black hole mass, declin- 
ing at higher black hole masses as the population is increas- 
ingly dominated by "dead" early-type systems which have 
decayed from their bright quasar epoch and may be fueled 
by quiescent (e.g. hot gas or stellar wind ac cretion) mecha- 
nisms. This trend agrees with that observed bv lHeckman et alJ 
(,2004), and predicts that black hole growth at z = is domi- 
nated by low mass systems, with Mbh 1O^M0. 

• The Mbh - o" Relation: Because feedback regulates ac- 
cretion and causes massive systems to enter a feedback- 
driven blowout phase in which the accretion rate declines 
and relatively little mass is gained, our model repro- 
duces and preserves an Mbh - cr relation similar to that 
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observed JGebhardt et all l2000t iFerrarese & MerrittI l200(i 
iTremaine et alJl2002^ ■ At the high mass end, where the ob- 
served Mbh - c relation is most robust, black hole masses are 
sufficiently large that any cold gas accretion described by our 
modeling will enter the "blowout" phase effectively immedi- 
ately, yielding no change in the Mbh - f relation. At lower 
luminosities, more detailed observations can test deviations 
from the high-MeH relation predicted by our modeling (see 
below), but we note that these are subtle effects and are con- 
sistent with all present observations of the relation down to 
quite low black hole masses. 

• Obscuration: The same "blowout" criterion which deter- 
mines the applicability of our feedback-driven solution and 
the Mbh - cr relation can be applied in a geometrical manner 
to determine the opening angle in which feedback will be suf- 
ficient to unbind gas and "punch through" unobscured sight- 
lines to the black hole. Objects with larger black hole masses 
will, in their initial high-m period, more efficiently unbind gas 
and blow out a larger fraction of gas, giving (indirectly) an 
anticorrelation between the obscured fraction and luminosity. 
In detail, our prediction of the obscured (Type 2) fraction as 
a function of observed (intrinsic) luminosity agrees with that 
estimated observationally from the same Seyfert luminosity 
functions we consider (Hao et al. 2005). 

Although qualitatively similar to "receding torus" mod- 
els w hich attempt to explain this effect (e.g., Lawrenc^ 
19911) or models of obscuration within a disk wind (e.g., 
Konigl & Karti dll994t ISto ne & Norrna ^994t iMurrav et alJ 
1^95; Elvis 200(1 iProga & Callman 20Q4|), our process op- 
erates over larger scales an d na turally explains the quan- 
titative distinction (see § 16. 1> between obscured frac- 
tions in local objects and (intrinsically) bright q uasars 
(e.g., lUeda et al."2003* 'Steffen et al.l"2003( iHasinged 1200411 
primes. RawHngs. & Willott 2004; Zakamska et al. 2004; 
ISimpsoni i2005^ Zakams ka et ab i2005i) . Our picture may, 
however, graft quite naturally onto such models (especially 
feedback-driven disk wind models) on smaller scales, as dis- 
cussed below. 

9.2. Testable Predictions 

There are a number of additional implications of our model 
for which observations either do not exist or are inconclusive, 
and therefore can be used to test our theory. 

• Evolution of Seyfert Luminosity Functions: From the evo- 
lution of the late-type galaxy luminosity function and gas 
fractions, we determine the evolution with redshift of the lu- 
minosity function determined by stochastic or steady-state 
cold gas accretion. Although there is substantial uncertainty 
in the late-type evolution, it is weak compared to the rapid 
evolution in the quasar luminosity function. Probing even 
moderate redshifts (e.g. z^O- 1-0.5) and separating the com- 
ponent of the quasar luminosity function contributed by late- 
type, non-interacting galaxies can test our prediction. We 
predict a relatively small and rapidly declining relative con- 
tribution to the faint-end of the quasar luminosity function at 
these redshifts, which has important implications for the con- 
tribution to cosmic backgrounds and buildup of black hole and 
spheroid mass, discussed below. 

• Eddington Ratio Distributions at Low Luminosity: Al- 
though we have compared our predicted Eddington ratio dis- 
tributions to observations as a function of mass, velocity 
dispersion, and luminosity, the observations are incomplete. 
Both the AGN luminosity function and Eddington ratio dis- 
tribution must turn over and begin to decline at some suffi- 



ciently low luminosity, and the large duty cycles calculated 
at the lowest Mb > -14 luminosities probed by e.g. lHao et alJ 
J2005h imply that observations are rapidly approaching this 
limit. Resolving the location and shape of this turnover, while 
nearly independent of the blast wave structure we have mod- 
eled (which determines the power-law slope of this distribu- 
tion above the turnover) determines the duty cycles and rate of 
excitation, constraining and testing our estimates of the con- 
tribution to luminosity functions from this fueling mechanism 
and its importance in the buildup of low-mass black holes. 

• The Mbh - c Relation: Because the mass of individual 
gas clouds is not correlated with a as is the total galaxy gas 
mass, and the timescale for "blowout" to occur is not deter- 
mined by the galaxy dynamical time in our scenario, we pre- 
dict a different slope for the Mbh - c relation than that indi- 
cated by models of spheroid and black hole co-formation (in 
e.g. g alaxy mergers) (e.g.. Silk & Rees 1998; Di Matteo et alJ 
120051) . At high black hole masses, black holes easily unbind a 
large gas mass rapidly and there is no change in the Mbh - o" 
relation, but at low masses, below an estimated break around 
^ 10^- IO^Mq, the slope of the Mbh-o" relation should be- 
come shallower, Mbh « for energetically-determined feed- 
back (or Mbh « for direct momentum feedback). This 
change in slope is a robust prediction, with a break set ap- 
proximately by the characteristic mass of a molecular cloud, 
but in detail the "break" location depends on uncertain fac- 
tors such as the efficiency of feedback coupling. Given this, 
it is possible that no change would be observable down to 
Mbh ~ IO^Mq, but even this observation would set strong 
lower limits to the feedback efficiency and means of feed- 
back coupling. For reasonable coupling values (similar to 
those which give the spheroid-formation induced relation), 
this gives a break at Mbh ~ lO^M©, and this slope change 
should be observable with larger samples of low-mas s black 
holes. In fact, existing observations by 'Barth e t all J2004 
2005) and Greene & Ho (2005a) favor a change in slope, but 
the large observational errors make this only a ^ 1 -2cr effect. 

From the observed spectrum of AGN column density distri- 
butions and molecular cloud sizes, we estimate their contribu- 
tion to the scatter in the Mbh - o" relation. Again, at high Mbh, 
they are negligible, but at low luminosities we predict that 
the intrinsic scatter in the Mbh - o" relation should increase 
substantially, from 0.27 dex to 1 dex at Mbh < lO^'M©. 
Present observations marginally la level) favor this, but 
expanded observations of the low-mass Mbh - o" relation can 
test this without any dependence on the relatively large sys- 
tematic uncertainties inherent in the absolute Mbh normaliza- 
tion as measured by different (e.g. reverberation-mapping vs. 
stellar dynamical) mass measurement techniques. 

• (Lack of) Global Effects on the Host Galaxy: A natural 
question following from our feedback-driven, blast wave sce- 
nario is whether or not there will be a substantial impact on 
the host galaxy. Indeed, for the case of gro wth-terminat ing 
"blowouts" in galaxy mergers (e.g., Hopkins et alJl20063) . a 
powerful wind expels or heats the remaining cold gas in the 
galaxy, enriching the X-ray halo, reducing column densities, 
and shutting down star formation. However, the situation is 
quite different in our modeling of cold gas accretion in quies- 
cent systems. In this case, black hole feedback is not forced 
to overcome the entire galaxy supply of cold gas, which is 
driven to the central regions and tightly bound by the entire 
galaxy potential via gravitational torques in a merger. Rather, 
only the relatively weak feedback to unbind cold gas clouds 
is imparted before blowout is entered, and most of the galaxy 
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ISM remains in an organized disk (which will not be impacted 
by the feedback). The feedback is essentially negligible as far 
as the cold disk of the host is concerned, although it may still 
enrich halo gas (as even a small quantity of metal-rich gas ex- 
pelled from the cloud "blowout" can be significant) and will 
heat diffuse bulge gas. There is no necessary causal connec- 
tion between this AGN activity and st ar formation, in t enta- 
tive agreement with observations (e.g..' Laine et al.l2006h . ex- 
cept insofar as larger quantities of cold gas will both increase 
the duty cycle for AGN excitation and enhance star forma- 
tion. The hosts will be relatively unperturbed, normal galaxies 
which continue to be actively star forming. 

• Testing Alternative Fueling Mechanisms: In addition to 
the predictions above, there are a number of tests which 
distinguish our model of accretion of cold gas and subse- 
quent feedback-driven light curve decay from other non- 
merger driven modes of accretion. For example, it has re- 
cently been suggested that stellar tidal disruptions can ac- 
coun t for the entire low-lumi nosity AGN luminosity func- 
tion ( Milosavlievic et al. 2006). However, stars cannot be dis- 
rupted outside the innermost stable circular orbit of the black 
hole for a black hole mass Mbh > 2 x lO'^M© (M^/Mq)"^ 
(where M, is the stellar mass). For the rates of stellar dis- 
ruptions estimated therein and by e.g.|Wan2 & Merritt (2004) 
to be sufficiently high to contribute the the observed lumi- 
nosity functions, the typical stellar mass disrupted must be 
< O.IMq, implying that the luminosity function is dominated 
by black holes with Mbh ^ 2 x XO^Mq (typically in Sd or 
Sm/Im galaxies). If more recent black hole m ass functions 
from Marconi et aP J2004l) : lSha nkar et al.' (^2^4) are adopted 
in estimating these rates, the expected black hole masses go 
down even further 

In contrast, our model predicts the luminosity function is 
dominated by ^ 10^ Mq black holes, in Sa/b and SO hosts. 
Furthermore, the bound mass in a disruption event is small, 
implying no correction to the slope or scatter of the Mbh - cr 
relation, and no association with cold gas supplies or out- 
flows (the entrained mass and momentum of which are much 
larger than can be accounted for by individual disruptions). 
Finally, the timescales for such events are very short - typ- 
ically ^ 1-lOyr. Although individual bursts of accretion 
may have shorter timescales than the time-averaged light 
curves we calculate, our characteristic timescales are consis- 
tent with both the integrated and "episodic" quasar Ufetime 
constraints from a number of observations (see e.g. iMartinil 
[2004, and references therein), which estimate lower limits to 
the episodic lifetime of > 10^ yr, and direct attempts to mea- 
sure the rates of short-duration X-ray flares (which hav e tra- 
ditionally been associated wit h stellar disruptions, e.g. iHiUsI 
119751 iMeszaros & Silklll977^ have found rates ~ 2-3 or- 
ders of magnitude below those needed to ac count for the 
low-lu minosity AGN luminosity function (e.g.. lDonlev et al.l 
l2()02l) . 

9.3. Comparison With The Bright Quasar Population 

Our modeling, in providing a galaxy-level context for the 
evolution of low-luminosity AGN in non-interacting galaxies, 
allows us to draw a physical distinction between the fueling 
of AGN and bright quasars, namely quiescent accretion vs. 
galaxy mergers and interactions. The feedback mechanisms 
and explosive "blowout" decay modes we determine here can 
be applied both to quiescent low-luminosity systems and to 
more violent, merger-driven systems (Hopkins et al- 2006 c), 
and on small scales (i.e. those relevant to detailed accretion 



and feedback mechanisms, molecular torii, and disk winds) 
the systems may be similar. However, there are fundamen- 
tal qualitative distinctions that should be kept in mind when 
making comparisons between the different classes of objects, 
and which imply quite different properties and physical inter- 
pretations of sometimes similar phenomena. 

• Fueling Mechanisms: The fueling mechanisms we 
have described here and in merger-driven activity as in e.g. 
[Hopkins et al. (2006a) are fundamentally different, the for- 
mer being determined by a steady cold gas supply in a qui- 
escent system, the latter by the violent torquing of cold gas 
throughout entire galaxies into the galaxy center in major 
mergers. Fueling by stellar winds or hot gas accretion may 
represent yet a third qualitatively distinct mode of fueling, 
which may not be be feedback-dominated but instead produce 
light curves dictated by steady-state equilibrium with hot gas 
(e.g., Ciotti & Ostrike^ [T997l" I2OOII) . or slowly decaying in- 
jection of gas from stellar winds inside the black hole radius 
of influence (e.g., Ouataert 2004; Soriaet al. 2005b). Fuel- 
ing by mergers is principally a cosmological quasar fueling 
mechanism, determined by galaxy-galaxy merger rates and 
the cosmological evolution in gas fractions and morphologies, 
characteristic masses of merging objects, and environmental 
dependencies of interactions. The fueling we have described 
here is non-cosmological, requiring only a local gas supply in 
disk galaxies, and independent of e.g. environment, formation 
mechanisms and times, and other cosmological influences (al- 
though indirect effects or correlations of these properties with 
e.g. morphology and gas fractions may be significant). 

• Luminosities: Objects fueled in these quiescent systems 
span a different range in luminosities from the bright quasars 
produced in gas-rich mergers. Our mechanism can account 
for local objects over the range -15 > Mb > -20, with an 
increasing contribution from post-merger decaying remnants 
above this, and mergers driving the brightest activity (with 
too low a number density to be seen at z = 0, but rapidly in- 
creasing to dominate populations down to Mb ^ -18 at mod- 
erate redshifts z > 1). Accretion of cold gas in quiescent sys- 
tems cannot explain the brightest quasars - in our derivation 
of the Mbh - c relation, it is clear that black holes cannot 
grow to the largest masses observed (compare the dashed and 
dot-dashed predictions of maximum masses grown via this 
method to the observed solid line at the highest black hole 
masses in Figure through this mode of accretion. Fur- 
thermore, even if this were possible, the high Eddington ra- 
tios at these masses needed to explain bright quasars cannot 
be maintained. Even with an extremely low coupling effi- 
ciency (further complicating attempts to match the Mbh - cr 
relation), a Mbh ^10^- 10^ M0 black hole would immedi- 
ately unbind the gas around it and yield a bright quasar life- 
time of ~ lO'*- 10^ yr, below observational lower limits (see 
e.g.. IMartin i 2004). The only way to force enough gas ac- 
cretion to "overpower" the blowout for a significant amount 
of time (enough time to have a chance to see the observed 
quasar population) would be to channel roughly the entire 
cold gas supply of the galaxy to the central regions at once 
- effectively requiring the tidal torquing possible only though 
a major merger. Major mergers account for high-luminosity 
activity naturally, with simulations producing typical bright 
quasar luminosities for the appropriate host galaxy masses 
(Hopkins et al. 2006a e). 

• Evolution: As emphasized above, fueling AGN via 
"stochastic" cold gas accretion as we model is fundamentally 
non-cosmological, whereas fueling via mergers is a cosmo- 
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logical process. This results in little evolution in the lumi- 
nosity function of objects driven by quiescent fueling, even 
assuming maximal evolution in disk galaxy populations and 
gas fractions up to z ~ 2. On the other hand, merger-driven 
activity evolves strongly with redshift, dominating above even 
a moderate Mb ^-19byz^ 1. The contribution from quies- 
cent fueling is quickly relegated to the faint-end of the lumi- 
nosity function at higher redshifts, becoming important only 
at luminosities one to two orders of magnitude 4 magni- 
tudes) below the break in the luminosity function at z 1, 
and another order of magnitude fainter relative to the break 
at z 2, the epoch of peak quasar activity. Thus, quiescent 
fueling is essentially irrelevant at high redshifts, while merger 
rates are high and merger-driven fueling dominates quasar ac- 
tivity. At higher redshifts, faint activity will be even more 
difficult to observe, but the distinction is also somewhat less 
clear, as essentially all systems are highly gas rich (and po- 
tentially unstable to gas collapse), but merger rates are also 
high. 

• Obscuration: Although the qualitative trend of a de- 
creasing obscured fraction with increasing luminosity is sim- 
ilar in both bright quasars and Seyferts fueled in the man- 
ner we m odel, there are important distinctions. As discussed 
in § 16.11 observations of local AGN and moderate-redshift 
populations (i.e. merger-driven quasar dominated, see Fig- 
ure O give a similar trend with luminosity but significantly 
different absolute obscured fractions. The Seyfert (local, 
low-luminosity) Type 2 fraction is systematically lower, ris- 
ing to comparable numbers of Type 1 and Type 2 objects at 
the lowest luminosities, whereas the quasar samples rapidly 
become dominated by Type 2 objects at the lowest lumi- 
nosities (^compare, e.g., | Hao et al. 2005; Uedaetal. 2003; 
[ Hasingedbool iGrimes. RawUngs. & Willott 2004: .Simpson. 
I2005h . Although both trends with luminosity are caused by 
black hole feedback, they are fundamentally distinct. In the 
local, non-interacting case the primary source of obscuration 
is small-scale, geometrical obscuration (see e.g.. lRisaiiti et^ 
[l999), and the feedback-based argument we use to describe 
the obscured fraction (giving e.g. comparable obscured and 
unobscured populations at low luminosity) is also a natu- 
ral consequence in many disk wind models and (potentially) 
torus models of small-scale obscuration. However, many X- 
ray identified bright quasars have been observed to b e grow- 
ing in starbursting or merg ing systems (e.g..; Alexand er et alJ 
l2005a'b:' Borvs et alJl2005 ). obscured not by a local geomet- 
rical structure but by the large-scale gas inflows powering ac- 
cretion, which can give a much larger obscured fraction at 
low luminosities and, in detailed simulations, gives a trend 
of obsc uration with luminosi t y whic h agrees with that ob- 
served ("Hopkins et al."2005b', '2006a'). This distinction be- 
tween quasi-static geometrical and strongly time-dependent, 
interaction-driven obscuration is important in understand- 
ing intrinsically bright Type 2 quasars (which are observed 
to be interacting, w ith galaxy-scale obscuring structures; 
IZakamska et alJl2005! ). and in building any model of obscu- 
ration for synthesis models of the cosmic X-ray or IR back- 
grounds. 

• Impact on the Host Galaxy: Merger-driven quasars are 
generally presumed to have a dramatic effect on their host 
galaxies, driving a galaxy scale-outflow, expelling gas and 
heatin g it to the virial temperature and terminatin g star forma- 
tion ( Springel et al. 2005a; Hopkins et al. 2006b), and enrich- 
ing the X-ray halo with metal-enriched gas (Cox et al. 20051. 
However, as discussed in §0 we do not expect the "blowout" 



from quiescent cold gas fueling to result in any significant im- 
pact on the host galaxy, even for small hosts. This can be un- 
derstood by considering the different thresholds for blowout. 
In the case of a merger, nearly the entire remaining gas mass 
is torqued to the central regions of the galaxy, in a roughly 
isotropic manner, and "blowout" is not entered until the en- 
ergy is sufficient to unbind this material. However, in the qui- 
escent case, the disk is largely unperturbed (and further has a 
limited covering angle), so it does not need to be "overcome" 
to enter the blowout - instead, only a small fraction of gas is 
unbound and it has little impact on the host. Recent observa- 
tions support this picture, with the outflow kinetic energy and 
entrained mass on large scales scaling with the galaxy size 
(and theref ore black hole mass) in a manner similar to what 
we predict (Bau m & McCarthv '2000). with Seyferts in small 
hosts doing little "damage" to their environment. This may 
still be significant for metal enrichment of the halo and con- 
tinued heating of the small amount of diffuse bulge gas, or 
e.g. the production of high-velocity clouds in the ejecta, but it 
does not disturb the normal, star-forming host disk. 

• Host Properties and the Mbh — cr Relation: As discussed 
above, systems fueled by cold gas accretion will be preferen- 
tially normal galaxies, SOs and Sa/b galaxies dominating the 
high-accretion rate population, whereas at higher luminosi- 
ties the population will mainly be (at the faint end) ellipti- 
cals decaying from previous merger activity and (at the bright 
end) objects still affected by major interactions. The M^n-a 
relationship as determined by these two fueling mechanisms 
may also differ, with the fact that the immediate gas supply 
in quiescent systems is not linked to the total gas supply of 
the galaxy (whereas it is in a strongly torqued merger), ulti- 
mately producing a somewhat shallower slope in the Mbh - cr 
relationship at low black hole masses (Mbh ^ 10^- 10^ M0), 
and increased scatter in the relationship at these masses. 

• Contribution to the Black Hole Mass Density and Cos- 
mic Backgrounds: Given the relative evolution of the lumi- 
nosity functions determined by mergers and interactions as 
opposed to quiescent or stochastic fueling, it is straightfor- 
ward to calculate the relative contribution to the black hole 
mass density and cosmic backgrounds. The rapid evolution 
in merger-driven quasar activity and its dominance at bright 
luminosities, especially during the period of peak quasar ac- 
tivity z > 1, implies that other fueling mechanisms will not 
contribute significantly to these quantities. Quiescent fueling 
in the manner we have modeled only adds significantly to typ- 
ical black hole masses at low mass (Mbh 10*' M©), as can 
be seen from comparison of the "increase" above the merger- 
driven Mbh- o" relation in Figure0 below the Mbh ^ 10^ Mq 
masses which dominate the black hole mass density. From en- 
ergetics, it is then clear that backgrounds such as the cosmic 
X-ray background, which can be entirely a ccounted for by 
merger-driven quasar activity (Hopkins et al. 2006a), receive 
little contribution from these alternative fueling mechanisms. 

9.4. Summary 

Our feedback-driven model for the evolution of accretion 
rates in systems which are non-interacting and fueled by 
cold gas explains a number of observations of AGN, and 
makes testable predictions for future observations. The mod- 
eling allows us to distinguish low-luminosity (stochastically 
cold gas accreting) "Seyferts" and high-luminosity (merger- 
driven) "quasars" in a physically meaningful manner based 
on their respective fueling mechanisms, even when both have 
similar high accretion rates. This distinction is critical in 
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understanding the properties of local objects, and their con- 
tribution (or lack thereof) to cosmological backgrounds and 
buildup of black hole and spheroid mass. It is also impor- 
tant to understand such distinctions when extrapolating prop- 
erties of local, quiescent objects to bright objects potentially 
fueled in mergers, which may be appropriate on small scales 
(e.g. disk properties and detailed accretion and feedback cou- 
pling mechanisms), but is inappropriate on larger scales (fuel- 
ing rates and mechanisms, obscuration mechanisms in at least 
some fraction of cases, and impact on the host galaxy). 

Our model for briefly fueling high accretion rate activity 
in quiescent, non-interacting galaxies is intended to provide 
a large-scale context, relating this fueling to the host galaxy 
and enabling a consideration of the statistics of such activity. 
We develop a self-consistent model for the accretion rate evo- 
lution within a wind or outflow, allowing the black hole light 
curve and accretion rate statistics to be calculated a priori. 
By no means does this model attempt to describe the small- 
scale gas thermal and ionization structure around the black 
hole. The local phenomena around the central engine, such 
as disk processes, the detailed coupling of feedback mech- 
anisms, ionization, broad and narrow-line region evolution, 
and pc-scale molecular torii, which have been (and continue 
to be) studie d in detail bo th theoretically and observation- 
ally (see e.g., Krolik 1999, and references therein), are ob- 
viously crucial to understanding AGN activity but are outside 
the scope of the large-scale fueling and feedback mechanisms. 
Ultimately, as long as feedback regulates AGN activity in a 
rapid, explosive manner, which we have shown for the case of 
cold gas accretion in quiescent systems requires only a small 
fraction (< 10"^) of the accretion energy to couple to the sur- 
rounding medium, the large scale evolution of fueling and ac- 
cretion rates should be described by our model. 

However, our work emphasizes the need to understand 
the connection between these small-scale phenomena, which 
are pre sumably universal in AGN unification models (e.g., 
lAntonu cci 1993), and large-scale models for fueling and feed- 
back in both quiescent systems as presented here or in in- 
teracting systems as in ^^opkins et al. (2006a). Although 
care must be taken in extrapolating properties of local qui- 
escent objects to interacting systems, the broad agreement 



between the predictions of feedback-driven blowout in both 
cases suggests that feedback coupling may be quite similar 
in both. Our model for feedback driven "blowout" and evo- 
lution of the Seyfert light curve grafts naturally onto disk- 
wind models in which the broad and narro w Une r egions are 
part o f a high-velocity ma gnetically (e.g . . iStone & NormanI 
fl99l . advectively (e.g.. lNaravan & Yil 11994. or radia- 
tively dShlosman. Vitello . & Shavivl Il985t iKonid & Kartid 
fT99l f Murrav et al. 1995t lElvisI 1200(1: IProgaetalJ 1200(1 
iProea & Kallman 2004) driven outflow from the accretion 
disk, which immediately suggest mechanisms for the cou- 
pling of feedback to the surrounding medium and a mapping 
from local scales (i.e. the accretion disk region from which 
the wind is driven) to large scales. Such a wind could di- 
rectly drive an outflow, with mass pile-up in a "snowplow" 
phase and the (initially radiatively, energetically coupled out- 
flow) powering a momentum-driven wind. Alternatively, it 
could collide with gas on cloud or central cold gas disk 
scales and shock, giving a pressure and energy source which 
drives a pressure-driven Sedov-type explosion. It is sugges- 
tive that these mechanisms manifest similarly, but there re- 
main many alternative possibilities for the coupling of feed- 
back. Studying this in the "laboratory" of local Seyferts, 
where AGN-d riven winds are well- known and studied in de- 
tail (see, e.g. 'Veille ux et al J 12005'. and references therein), 
may potentially represent the ultimate mapping between the 
well-studied properties of AGN on small scales and the large- 
scale implications for the host galaxy and AGN statistics as 
studied here. 
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